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1.  Introduction 

In  this  paper  we  investigate  the  stability  of  equilibrium 
of  several  homogeneous  isotropic  elastic  solids.   The  problems 
considered  are  treated  by  making  use  of  a  non-linear  three- 
dimensional  theory  of  elasticity  of  the  same  sort  as  treated  in 
the  book  of  Green  and  Zerna  [l],  although  the  formulation  of  the 
theory  used  here  follows  that  of  Murnaghan  [2].   No  restrictive 
assumptions  are  made  concerning  the  thickness  of  the  bodies,  and 
the  displacements  and  strains  are  not  assumed  to  be  small.  The 
strain  energy  function,  from  which  the  stress-strain  relations 
come,  could  in  principle  be  chosen  arbitrarily;  however,  in  this 
paper  a  special,  but  for  many  materials  reasonable,  choice  for 
the  strain  energy  function  is  made. 

In  each  problem  a  simple  solution  is  derived  which 
satisfies  all  conditions  of  the  exact  non-linear  theory.  A 
small  perturbation  with  respect  to  this  simple  state  is  then 
made,  and  the  question  is  raised  as  to  whether  other  possible 
equilibrium  states  —  buckled  states  in  other  words  —  exist 
within  the  framework  of  the  basic  theory.   Critical  pressures 
and  strains  needed  to  cause  such  buckled  solutions  are  found, 
together  with  the  corresponding  modes  of  buckling.   It  is 
observed  that  the  critical  buckling  pressures  and  strains  are 
the  same  as  those  furnished  by  the  classical  thin  body  theory 
in  the  limit  case  of  small  thickness. 

Below  we  list  a  brief  description  of  the  problems 
considered. 


Problem  A:   A  circular  cylinder  is  compressed  along  its 
curved  lateral  surface.   The  compression  is  such  that  the  curved 
lateral  surface  of  the  undeformed  cylinder  goes  into  the  curved 
lateral  surface  of  a  coaxial  circular  cylinder  of  smaller  radius 
in  such  a  way  that  no  shear  stress  is  developed.   It  is  as  though 
the  cylinder  were  compressed  by  shrinking  a  very  stiff  greased 
ring  down  on  it.   The  plane  faces  are  assumed  to  be  free  of 
stress.   We  look  for  buckling  in  the  axial  direction. 

Problem  B;   A  hollow  circular  cylinder  has  a  hydrostatic 
pressure  applied  to  its  outer  curved  lateral  surface.   The  inner 
curved  lateral  surface  is  assumed  to  be  free  of  stress.  The 
axial  displacement  is  taken  to  be  zero.   Thic  can  be  brought 
about  by  requiring  that  the  plane  faces  of  the  cylinder  lie  on 
rigid  greased  plane  plates,  the  distance  between  the  greased 
plates  being  the  unstrained  height  of  the  cylinder.   VJe  look 
for  buckling  in  the  plane  of  the  faces  of  the  cylinder. 

Problem  C:   A  circular  cylinder  is  subjected  to  axial 
compression.   The  compression  is  such  that  the  end  faces  remain 
plane  and  horizontal  and  no  shear  stress  is  developed.   This 
might  be  brought  about  by  compressing  the  ends  of  the  cylinder 
between  tv;o  rigid  plane  horizontal  greased  plates.   The  curved 
lateral  surface  of  the  cylinder  is  assumed  to  be  free  of  stress. 
We  look  for  buckling  perpendicular  to  the  axis  of  the  cylinder. 

Problem  D:   This  problem  is  the  same  as  Problem  C  except 
that  the  cylinder  is  hollow  and  the  inside  surface  is  assumed 
to  be  free  of  stress. 


Problem  E:  A  hollow  sphere  has  a  hydrostatic  pressure 
applied  to  its  outer  surface.   Its  inner  surface  is  assumed  to 
be  free  of  stress. 

Problems  B  and  E  have  been  treated  previously  by  Lubkin 
[3]  using  somewhat  different  methods.   The  other  problems  have 
also  been  considered  previously,  but  we  believe  not  with  the 
general  non-linear  theory  used  here. 

2.  The  Non-Linear  Theory 

In  this  section  we  outline  the  theory  to  be  used,  largely 
without  derivation.   The  notation  used  is  that  of  Fritz  John 

W,   [5]. 

Consider  a  rectangular  Cartesian  reference  frame  X.  A 

particle  which  is  at  the  point  (x,,X2,x^)  when  the  body  is 

unstrained  will  be  at  a  point  (x,,Xp,x^)  after  the  body  is 

strained.  The  Eulerian  coordinates  x, ,  Xp,  and  x,  are  treated 

as  functions  of  the  Lagrange  coordinates  x-,  ,  Xp,  and  x^. 

Dependence  on  time  is  ignored  since  we  consider  only  equilibrium 

problems.       _ 

^x. 

Let  p.  .  =  -c-r^  (i,j  =  1,2,3).   Then  W,  the  strain  energy 
ij   ox. 

per  unit  undeformed  volume,  is  a  function  of  the  quantities  P^  .:  • 

Let  q.  ..  =  ^^      (i,j  =  1,2,3).   If  v;e  neglect  body  forces, 
the  equations  of  equilibrium  are  [5] 


■3r 


(2.1)  -^  =  0  ,  i  =  1,2,3 


where  we  use  the  usual  summation  convention. 


Let  S  and  S  denote  surfaces  containing  the  same  particles 

in  the  unstrained  and  strained  bodies  respectively.   Let 

-»  ±         —    —     — 

n  =  (n, ^np,n^)  and  n  =  {n,,n2,n^)  be  unit  normal  vectors  pointing 

to  corresponding  sides  of  S  and  S.   Let  t  =  (t-|^,t2,t,)  be  the 

surface  traction  vector  exerted  on  the  surface  S  from  the  side 

to  which  n  points  (the  units  of  t^  are  force  per  unit  deformed 

area ) .   Then 

(2.2)  t^dS  =  qij-n^dS  ,  i  =  1,2,3 

where  d'S  and  dS  are  the  elements  of  area  of  S  and  S  respectively. 

We  can  think  of  the  transformation  dx^  =  p^  .dx .  as 
consisting  of  a  rotation  and  a  pure  deformation  which  are  defined 


(^^p)""^ 


as  follows.   Let  p  be  the  matrix  (p^.  .),  and  c  =  p(/p*p 
(*  denotes  the  transpose  and  ^p^p  is  the  symmetric  positive 
definite  square  root  matrix  of  p*p).   Then  c  is  orthogonal  and 
det  c  >  0  (we  assume  det  p  ^  0)  so  that  c  is  a  rotation  matrix. 
For  reasons  stated  later  we  call  c  the  local  rotation  matrix  of 
the  deformation  [6].   We  define  pure  deformation  matrices  A  and 
B  by  the  equations  p  =  Ac  =  cB  (i.e.   p  is  the  rotation  followed 
by  a  pure  deformation  and  p  is  a  pure  deformation  followed  by 
the  rotation).   We  continue  with  a  discussion  of  the  rotation 
and  pure  deformation  matrices. 

Let  ds  and  d'S  be  the  differentials  of  arc  length  along  a 
curve  consisting  of  the  same  particles  in  the  undeformed  and 

deformed  bodies  respectively.   Then  the  eigenvalues  of  p*p  and 

d'S  2 
pp*  are  the  stationary  values  of  (^)   with  respect  to  variations 

of  direction.   The  eigenvectors  of  p*p  and  pp*  give  the 


ds 
directions  of  stationary  ^  in  the  undeformed  and  deformed 

bodies  respectively.   The  matrix  p  carries  each  eigenvector  of 

p*p  into  an  eigenvector  of  pp*.   Thus  if  z  is  a  unit  eigenvector 

d's 
of  p*p  and  A  is  the  corresponding  stationary  value  of  -r— ,  then 

Q.O 

z"  =  Y  P^  is  a  unit  eigenvector  of  pp*.   Since 

cz  =  p(yp*p)   z  =  T-  pz  =  z",  we  see  that  c  rotates  the  dlrections- 

of  stationary  ^  in  the  undeformed  body  into  the  directions  of 

d's 
stationary  -r—   in  the  deformed  body.   This  is  the  reason  we  call 
ds 

c  the  local  rotation  matrix  of  the  local  deformation  p. 

Since  c*p  =  /P*p,  we  see  that  c*p  is  symmetric  and 

positive  definite.   We  next  prove  that  these  conditions 

characterize  c.   Let  d  be  any  other  rotation  matrix  such  that 

d*p  is  symmetric  and  positive  definite.   Then 

(d*p)   =  (d*p)*d*'p  ==  p*p  so  that  d*p  is  a  square  root  matrix 

of  p*p.   But  there  is  only  one  symmetric  positive  definite 

square  root  matrix  of  p*p,  namely  c*p.   Hence  d*p  =  c*p  and 

d  =  c. 

-1 


Similarly  C  =  (/pp*)   p  is  the  unique  rotation  matrix  such 
that  pC*  is  symmetric  and  positive  definite.   However 
Cz  =  (  ypp*)   pz  =  A(7pp*)   1  =  "z   for  all  unit  eigenvectors  z 
of  p*p.   Hence  C  =  c  and  c  is  also  characterized  by  the 
condition  that  c  is  a  rotation  matrix  and  pc*  is  symmetric  and 
positive  definite.   Also  pc*  =  Tpp** 

From  the  previous  definitions  of  the  pure  deformation 


matrices  A  and  B  we  see  that  A  =  pc*  =  Jpp*  and  B  =  c  p  =  yp*p. 
Hence  the  eigenvalues  of  the  pure  deformation  matrices  are  the 


ds 
stationary  values  of  -r-,    and  the  eigenvectors  give  the  directions 

Q.S 

ds" 
of  stationary  ^  in  the  deformed  and  undeformed  bodies 

respectively. 

Natural  candidates  for  the  strain  matrix  are  the  pure 
deformation  matrices  minus  the  identity.   We  define  r)  =  J  p*p  -l 
to  be  the  strain  matrix. 

We  now  list 

(2.3)    cc*  =  1,  det  c  >  0,  and  c*p  (or  equlvalently  pc*) 
is  symmetric  and  positive  definite. 


(2.4)  T)  =  Jp*p-  1  =  c*p  -1  . 

The  matrix  c  denotes  a  rotation  around  an  axis  through  an 
angle  which  we  call  the  local  rotation  angle.   Let  (r, ,rp,r,) 
be  a  unit  vector  along  the  axis  of  rotation  and  denote  the  local 
rotation  angle  by  ij/.      f   is  taken  to  be  positive  if  the  rotation 
c  would  cause  a  right  handed  screw  to  advance  in  the  direction 
of  (r-|^,r2,r^).   Then 

(2.5)  c..  =  5..cos^  +  r.r.  (l-cos^)-e..,  r,  sin^ 

where  e.  .,  =  ±1  if  i,  j,  k  is  an  even  or  odd  permutation  of 
1,  2,  3  and  e.  .,  =  0  if  two  subscripts  are  equal.   In  some 
problems  considered  here  one  can  determine  the  quantities  r.  by 
inspection  and  then  obtain  ^  as  a  function  of  the  quantities  p 


10 


using  the  symmetry  of  the  matrix  c*p. 

For  an  isotropic  material  the  strain  energy  density 

function  W  should  be  taken  as  a  symmetric  function  of  the 

d^ 
stationary  values  of  -3—.   Hence  for  such  a  material  W  is  a 

as 


function  of  s-,  ,  Sp,  s,  where  s.  is  the  sum  of  the  i-th  powers  of 
the  eigenvalues  of  ri.   If  we  let  the  square  bracket  denote  the 
trace  of  a  matrix,  then 

(2.6)  s^  =  [nM  =  [(c*p-l)M  . 

To  obtain  an  expression  for  the  Lagrange  stresses  q. .  we 

as,  ^' 

derive  expressions  for  the  partial  derivatives  ■^- — .   Since 


?!f^* 


3s,  J" 

ij    ""    "IJ 


s^  =  M   =  [c*p]  -  3,  we  have  -^  =  c^.  +  [-^^   p].   Prom  (2.3) 
we  obtain  ^^^ —  =  -c*  -5- —  c*.  Hence 

tt^  p]  =  -t^*  ^  ^*p^  =  -f^^*  ^  ^^)*1  =  -^^'^  tr-  ^^ 
..nf^c)(p*c)]  =  -[||^cc*p]  =  -[|£^p]  . 

Thus  tlfrrP^  =  0  ^"^-357-^  ^ij- 

Next 

Sg  =  [n^]  =  [t*^]  =  [(p*c-l)(c*p-l)] 

=  [p*p]  -2[c*p]  +  3  =  [p*p]  -2s^  -  3 

so  that  ^=  2(p   -c   ). 
Finally 
s,  =  [n^]  =  [(P*c  -  l)(p*p  -  2c*p  +1)] 

=  [p*cp*p]  -  3[p*p] +  3[c*p]  -  3  =  [p*cp*p]-  3(S2  +  S3^  +1)  , 
and 

-4^  =  ^^   [p*cp*p]-  6p   +3c  .  . 
op.j   op.j  ij    ij 
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The  procedure  for  differentiating  [p*cp*p]  is  similar  to  that 

for  [c*p]  and  the  result  is  -5- —  [p*cp*p]  =  3p.-,c  p   .   Thus 

t'P.  .  ir  SI  sj 

OS  ?^ 

BpTJ  =  ^Pi^c^^p^  .  -  6p.  .  +  3c,  ..   Since  q^^  =  ||-  ^  we  have 

therefore: 

q.  .  =  (  -"v 2  -5 1-  3  -T —  I  c .  . 

^ij   V°^i    °^    osZ  y  ij 

(2.7) 

+  H"Si;-^^>ij''^^Pir^srPso  • 

This  derivation  of  the  quantities  q. .  was  shown  to  the 
author  by  F,  John;  although  this  author  is  not  aware  that  it 
appears  in  any  of  P.  John's  writings. 

The  special  strain  energy  function  used  here  from  now  on 
is  what  might  be  called  the  standard  strain  energy  function: 

(2.8)  ^  "  ■?  ^1  ■^'-''^2  * 

It  corresponds  to  the  classical  strain  energy  density  function 
when  A  and  \i   are  the  Lame  constants.   For  this  strain  energy 
function,  (2.7)  becomes 

(2.7')         q^.  =  (Xs^-2^i)c..+2,xp,. 

and  the  equilibrium  equations  for  zero  body  force  become 

Sp.  .  be.  .  Ss, 

(2.9)  2^.^+(As,-2^.)  ^+Ac^^=  0  . 


lZi\0  . 


3.  Introduction  of  Tensor  Methods 

In  the  problems  treated  here  we  find  It  convenient  to 
introduce  tensor  methods  and  work  directly  in  curvilinear 
coordinates. 

Let  9     =  9.(x,,X2,x,)  (1  =  1,2,3)  be  the  curvilinear 

bx     hx 
coordinates,  and  let  g.  .  =  -^g^  -^^   be  the  9-components  of  the 

1   J 

metric  tensor.   The  g"''"^'s  are  defined  by  g  g  .j^  =  5,  ,  and  the 

g.  .'s  and  g  "^ '  s  are  used  in  the  usual  manner  to  lower  and  raise 

indices  of  the  ©-components  of  tensors. 

In  Table  1  the  right  hand  column  presents  symbols 

representing  the  0-components  of  tensors  whose  X-components  are 

represented  by  the  symbols  in  the  left  hand  column.  For  example, 

^x  Sx 
^11  ~  "5^  "3^^  ^rs*   ^®  call  special  attention  to  the  fact  that 

the  tensors  considered  have  the  items  in  the  left  hand  column 

as  X-components  rather  than  as  X-components. 


X-component£ 

3 

0. 

■components 

^i 

"i 

PiJ 

^j 

^Ij 

«ij 

■^i 

"i 

^1 

\ 

h 

\ 

^0 

"u 

^1 

\ 

®ijk 

Table 

1 

^IJk 
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To  give  physical  meaning  to  some  of  the  quantities  intro- 
duced by  Table  1,  we  observe  the  following.   Consider  the 
vectors  x  =  (x,,Xp,x,),  x  =  (x,,  x^,  x.,),  n,  n,  and  t.   Let 

g.  =  -T^-  so  that  g.  is  a  tangent  vector  to  the  i-th  curvilinear 
curve .  Then  x  =  vx'^'g,  ,  n  =  N  g .  ,  n  =  N^g".  ,  and  t  <=  T  "g .  .  Thus 
u  ,  etc.,  are  the  components  of  the  corresponding  vectors  with 
respect  to  the  g^'s. 

Also  we  observe  that  P.  .  =  u.   where  I   denotes  covariant 


ij   "i 


0       'J 


differentiation  with  respect  to  9  .  and  e.  .,  =  /g  e.  .,,  v;here 
g  =  det  (g^j)' 

In  curvilinear  coordinates  some  of  the  results  of  Section 
2  become 


(3.1)   Q^J|j  =  0  , 


ij 
(5.2)   T^dS  =  Q^^N  dS  , 

J 

(3-5)   Cj^^C^-^  =  g^"^,  det  (Cj,  )  >  0,   (C^^Pj^*^)  is  symmetric  and 
positive  definite, 

(3-5)   C^  =  g^j  cos  ^  +R^Rj(l  -  cos  f)  -  Eiji^R*^  sin  ^  , 

(3.6)   s^   =  cJ^Pj^  -3,  etc.. 


^^'^'    ^ij  A^"^^''^^;^! 


3 

+  2r|^-3l^V,,  +  3|^P,  c^^P  . , 

V^Sg    osZJ    ij   ^  ds,  ir    sj 


(3.7')  %■   =   (As,  -  2ij,)C   +  2nP.  .  for  standard  strain  energy  , 


11 


. .  . .  . .  s?, 

(3.9)   2M.P  -^1  +  (As^  -  2|i)C^J  I  _  +^C^'^  ^  =  0  for  standard 
strain  energy. 


4.  The  Perturbed  Problem 

In  each  of  the  problems  considered  we  assume  that  the 
buckled  solution  can  be  expressed  as  a  function  of  a  parameter 
5^0  such  that  the  buckled  solution  coincides  with  the  simple 
solution  for  5=0  but  differs  from  the  simple  solution  for 
5  >  0.   We  also  assume  5  can  be  chosen  so  that  quantities 
depending  on  6  have  a  derivative  with  respect  to  6  at  5  =  0. 
If  A  is  any  such  quantity,  we  let  A  be  the  value  of  A  for  6=0 
and  A  the  derivative  with  respect  to  5  at  5  =  0. 

In  the  perturbed  problem  we  treat  the  u  ' s  as  the  unknown 
functions.   The  main  object  of  this  section  is  to  derive 
differential  equations  for  the  u  's. 


First  of  all  from  P. .  =  u. I   we  obtain 

1 J    11  i 

(4.1)  P,  ,  =  u, 


3 
Prom  (3. 3)  we  have 

('4.2)  c^^b^^  +b^^cj   =  0  , 

('+.3)         {t^HJ^+C^^?J)   is  symmetric  . 

Since  the  left  side  of  equations  (4.2)  are  the  elements  of 
a  symmetric  matrix,  (4.2)  can  be  thought  of  as  six  linear  equa- 
tions for  the  C'^'^'s.   Condition  (4.3)  gives  us  three  more  linear 


■„>n  ..  «•-■      .  j  : 


■  ^     ,.  '--^•ir|p  ""J 


■'i : 
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'11 
equations  for  the  C  '-''s,  so  (4.2)  and  (4.3)  together  determine 

the  C'^'^'s  as  functions  of  the  P'^'^'s  and  hence  as  functions  of  the 

u^'s  through  (4.1). 

Next  we  derive 


(4.4) 


s. 


1^2 


«     •  1  T 


S^  = 


iki. 


^(^iJ-^Pij-^Pi^C-^^kj^^''^ 


;,i. 


which  give  the  s. 's  in  terms  of  the  u  's. 

Since  s  =  [c*p]-3,  we  have  s^  =  [c*p]  +[c*p].  We  can  show 


that  [c  p]  =  0  in  the  same  way  we  showed 

Section  2.   Hence  s^  =  c^^^  .   =  C^.P  "^ . 
Since   Sp  =  [p*p]  -2s,  -3^  we  have 


^ 


[c*p]   =  0  in 


ij 


Since      s^  =    [p*cp*p]  -  3( Sp  + s^  + 1 ) ,   we  have 

s^  =   [p*c  p*^]  +  3[pc    pp   ]  -  3(2p^j  -c^j)p^j    . 

Methods  similar  to  those  used  to  show  [c*p]  =  0  can  be  used  to 
Show  that  [p*cp*p]  =  0.  Thus  s^  is  given  by  (4.4). 

By  dotting  (3.1)  and  (3.7)  we  now  obtain  the  differential 
equations  for  the  u  's  and  the  Q  '^'s.   Since  we  are  using  the 
standard  strain  energy  function  in  our  problems,  we  list  these 
results  only  for  that  special  case. 
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(4.5) 
(4.6) 


Q^"^  =  (aI^-  2m.)C^'^"  + As^C^"^'  +2m.P^"^  , 


aj 


a's- 


Ss. 


=  0 


5.  Special  Curvilinear  Coordinates 

In  the  problems  treated  here  either  cylindrical  or 
spherical  coordinates  are  used.   In  this  section  we  list  some 
results  for  these  two  curvilinear  coordinate  systems. 

For  cylindrical  coordinates  we  have  x,  =  0,  cos  0 
Xp  =  d     sin  0       x^  =  ^^'      For  convenience  we  write  r  =  9,  , 
9  =   9^,    z  =   9         Then 


X. 


V  Xn  ,  Xp  J  X^  / 


7 "^2 


/ 


/ 


Fig.  1 


1^ 


(5.1) 


r 


{ 


1 

f  1 

0 

0 

'%j'  =( 

0 
V  0 

0 

0 

1 

(s^J)  '[ 

^1 

0 

0 

0 
0 

\ 

\0 

0 

1 

=  r 


The  non-zero  Chrlstoffel  symbols  are 


(5.2) 


^2  2 


12 


21   "   r 


r 


22 


=  -r 


Since 


and 


we  have 


'\  - "' 


,iJ 


,i 


^  +    '    kj^        ■"   '    k/        ' 


hu^        2     au' 


(5.3)  (P^j)=|     if^u^       |#4ul       ^' 


Bu^ 


\'Sr 


^5~ 


<3z 


r 


(5.4) 


,1J 


>2j 


^^^1.= 


p3J| 


j 

5?^-^'    ,    1  p31 
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For  spherical  coordinates  we  have  x,    =   9     sin  9     cos  9 
Xg  =  Q-j.   ^^^  ^2   ^^^  ^5'    ^'5  ~   ^1   ^'^^   ^2*      ^°^  convenience  we  write 


=   9^,    9  =  9^,    ^  =  9 


\  X-,  fX^fX-z) 


>x. 


Fig.    2 


r 


(g. .)  = 


<    (g'J")    = 


10  0 

0  r^       0 

0         0       r^  sin^e 


/ 


/l 
0 

\  0 


1 

0 


0 
0 


r      sm   9 


g   = 


r      sin   9 


The  non-zero  Christoffel  symbols  are 


r"  ji^l 


(5.6) 


r 


22 

,2 
12 


=  -r 


,      r^^  =  -r  sin  0  ,   P 


33 


33 


=  -sin  9   cos  0  , 


1  21    I  1- 


=  r 


31 


r 


r 


23 


r-.3 

I  ^2  =  cot  0  , 


also 


16 


^^  -ru^sm^e 


(5.7)  (pV  = 


du^    1    2    au^  ,  1  ,1 


-u^sln  9  cos  9 


or   r     60 


cot  0  Ix-"^-?  u-'-+u^  cot  9, 


V^"l 


(5.8)  /  p 


2j 


oPii  ^ 2  pll _  ^p22  _  ^p55  3i^2  , ^  pl2  ^^^  g 

69 .    r 
J 

Sp2J 


~  "^^ 


+  P^^  cot  9  -  p55  sin  e  cos  9  +  i  (P^^+^P^^) 


.3j 


|P_i  +  i  (pl5+ 3P^1)+  (p25  +  2P^^)  cot  9 

J 


6.  Problem  A 

This  problem  was  treated  earlier  by  this  writer  [7]  as  a 
separate  problem.  It  seems  advisable  to  include  it  here  since 
this  treatment  is  considerably  simpler. 

Cylindrical  coordinates  r,  9,  z  are  used.   The  unstrained 
circular  cylinder  occupies  the  region  r  _^  R,  Izj  _^  h. 


^  X, 


Fig.  3 
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We  observe  that  ^,=  (cos  9, sin  0,0),  gp =  r( -sin  0,cos  9,0) 
and  g-z  =  (0,0,1)  (see  Section  3).   Since  gg  and  g^  are  per- 
pendicular to  the  radial  direction,  the  condition  that  the 
curved  lateral  boundary  goes  into  the  curved  lateral  boundary 
of  a  coaxial  cylinder  with  smaller  radius  may  be  written  as 
u  =  (l-e)R  for  r  =  R  where  e  is  a  positive  constant.   Since 
gp  and  g,  are  tangent  to  the  deformed  curved  lateral  surface, 

the  condition  that  no  shear  stress  is  developed  for  r  =  R 

2    3 
becomes  T  =  T^  =  0.   Hence  the  boundary  conditions  can  be 


written  as 

for  r  =  R 


T""  =  T^  =  0  j 


(6.1)       <  ^"''  =  (l-e)R  j 


T^  =  0   (i=  1,2,3)   for  z  =  ±h 


v^ 


On  the  boundary  r  =  R,  N,  =  1,  N^  =  N,  =  0  and  on  z  =  ±h, 
N,  =  Np  =  0,  N^  =  ±1.   From  (3.2)  we  see  that  the  boundary 
conditions  can  be  rewritten  as 


Tq^I  =  Q^l  = 


(6.2) 


We  look  for  a  simple  solution  in  which  the  cylinder  is 
linearly  extended  in  the  radial  and  axial  directions  with  no 

rotation  around  the  axis  of  the  cylinder.   Since  x  =  u  g^,  this 

"I         2        "^ 
means  u  =  ar,  u  =  0,  u-"^  =  bz  for  the  simple  solution  where  a 

and  b  are  constants. 

From  (5.3)  we  have 
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Since  (P  '^ )  is  symmetric  and  positive  definite  (for 
a,b  >  0),  we  have  C^^  =  g^J  from  (3.3). 

From  (3.6)  we  have  s,  =  2a+b-3. 

Substituting  P-'""^,  C"^-^,  and  s^  into  (3.9),  we  see  that  the 
simple  solution  satisfies  the  equilibrium  equations. 

From  (3.7')  we  see  that  (Q"'"'^)  is  diagonal. 


and 


Q-^-"  =  2(  A  +  n  )a  -  3A  -  2u  ,   Q   =  -^  Q   , 

r 

Q^^  =  2Aa  +  ( A  +  2p.  )b  -  3^  -  2(1  . 
Substituting  into  (6.2)  we  obtain 


I  a  =  1-e 
(6.3)  I 

■  b  =  1  +  .  ""X  =  1  + 


I  b  =  1  +  -,  ,  ^  =  1  +  •= 


where  (T  is  Poisson's  ratio. 

Thus  we  observe  there  is  a  simple  solution  of  the  type 
considered  v;hen  0  j^  e  <  1. 

On  the  boundary  r  =  R,  d S  =  abr  drd9  and  dS  =  r  drdQ. 
Hence  letting  T  be  the  normal  pressure  per  unit  deformed  area 
applied  to  the  edge  of  the  plate,  we  obtain 

T  =  T^  =  o^^  1^  -  gj^  (  3A  -t-  2u )  £  _      Ee 

~    dS     A  +  2|i    ab   (1-e)  ( 1  - -T+  2re  ) 

where  E  is  Young's  modulos.   Hence  T  — >  00  as  the  radius  of  the 
deformed  cylinder  goes  to  zero  (i.e.  as  e  — >  1). 
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If  we  let  P  be  the  total  load  applied  to  the  edge  of  the 
deformed  cylinder,  we  have  P  =  4-jrabRhT  =  ^-jrRhe  -^ — -  .   Hence  the 
total  load  approaches  a  finite  limit  as  the  deformed  radius  goes 
to  zero.  Also,  as  the  radius  goes  to  zero,  the  deformed  height 
2bh  approaches  a  finite  limit.   Hence  a  finite  total  load 
shrinks  the  cylinder  down  to  a  line  segment  of  finite  length. 
From  these  considerations  we  see  that  the  standard  strain  energy 
function  does  not  give  a  good  description  of  known  materials 
unless  the  strains  are  not  too  large.   However,  no  limitations 
on  the  magnitudes  of  strains  are  assumed  in  the  following. 

We  look  for  special  buckled  solutions  which  are  axially 
symmetric,  and  therefore  are  free  of  displacement  in  the  0- 
direction.  Hence  vie   take  u  =u(r,z,6),u  =0,  u-^=  u-^(r,z,5) . 
We  write  u  =  ar +V  ,  u^  =   bz+V-^.  From  (6.3)  a  and  b  are 

functions  of  e  which  in  turn  is  a  function  of  the  perturbation 

•  1     *  ^ 
parameter  6.  Hence  u  and  u  (see  Section  4  for  the  meaning 

of  "^  and  '  )  have  the  form 

00 


u-^  =  ar  +  >   a„(z)a  (r) 

(6.4)       { 

00 

u^  =  bz  +  ^~~;  b„(z)0  (r) 
^ — ?T  n    n 
n=0 

"1     *  "5 
where  we  have  used  series  for  V  and  V-^.  That  is,  we  express 

•  1     •  ^ 

u  and  u-^  as  the  perturbed  simple  solution  plus  a  sum  of 

products.   We  will  attempt  to  choose  the  «  's  and  P„'s  so  that 

variables  separate  when  u  and  u^  are  substituted  into  (4.6). 


From  (5.3),  (4.2),  (4.3),  and  (4.4)  we  have  (remember 

O  4        0        J,  A 

J   ^ij'   j    J    J 


C    -  g   ,  C    -  g^  ,  C   «C   _  b  } 
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(6.5) 


r 


(pij)  = 


0 


dz   \ 


au5 


2l 


a  +b 


t   523.  _1_  (p23.p2j  ^  Q 
a  +b 

s-|^  =  2a  +b  -  5 


Su^  .  1  -1  .  Su^ 


^1  =  ^F--^F  ^  +^ 


From  (5.^),  (6.5),  and  (6.^) 


and 


plj 


p3J 


[a  (a"+i  a'  -^  a  )  +  a"a  ]  , 
■■  n  n  r  n  „2  n    n  n-*  ' 


1  o» 


P^'^l  =  r~  [h   (P"  +i  P')  +b"p  ] 
Ij   ^ n  n  r  n    n'^n 


To  separate  the  variables  in  (^.6)  we  try  to  choose  a  and 


P  so  that  the  variables  separate  in  P  "^  |  .  and  P  '^  I   themselves, 

■^         J 

I.e.  we  choose  a.     and  P  so  that  a"  +—  a'  — —  a     and  a  are 

n      n  n     T     n       d     n  n 

T  r 

proportional  and  so  that  P!! +-7:  P'  and  P  are  proportional.  This 
can  be  accomplished  by  choosing 


(6.6) 


a^   =  Jl(  V^ 


L^^'^o^v) 
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where  J  and  J,  are  the  usual  Bessel  functions,  k  =  0  and 
o      1  '  o 

k  "^   ^n+l  ^°^  ^  ~   0>1>2,...  are  constants  to  be  chosen  at  our 
convenience. 

Since  k  =  0,  we  see  that  a  =  0  and  P  =  1.  It  follows 
from  (6.4)  that  a  Is  arbitrary.  For  convenience  we  choose 


^o  =  °- 


We  now  have  the  following  relations  for  a  and  P  : 


(6.7) 


^^0"+!  a'  +  (kj-i)a  =  0  , 


n  r  n 


n  „c:  n 


P"  +^  P'  +k^p  =  0  , 
n  r  "^n   n  n     ' 


n  r  n 


P  =  -k  a 
•^n     n  n 


n  n 


The  first  two  equations  are  the  differential  equations  for  J, 

and  J  and  the  last  two  are  familiar  recurrence  formulas  for  J. 
o  c 

and  J, .  These  can  be  used  to  replace  derivatives  of  a  and  B 
by  expressions  involving  a  and  P  themselves. 


Let  c^  =  a;+k^b^  and  d^  =  b^^  Vn" 


Using  (6.5),  (6.4) 


and  (6.7)  we  have 


and 


0^5  =  -i-  ri  c  a 


S-,  =  2a  +  b  + 


d  P 
n  n 


In  the  same  way  it  follows  that 
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I  =  T~Z   (c'  -k  d  )a    , 

I  i   ^ n   n  n  n   ' 

Ph  =  V—   (d'  -  k  c  )(3 

Ij   ' n   n  n  n  * 

'  J   2  .  £  "^^ n  n  ' 


a  +b 

'o     ^+b  ^ "  '^  ^ 

Next  observe  that  aI^  -  2m,  =  -2y,b  and  substitute  the  above 
into  {k.6).     The  variables  separate,  and  we  obtain  the  following 
ordinary  differential  equations. 

[Bc^-  (A+2n)k^d^  =  0  i 
(6.8)  ^  V   n  =  0,1,2,... 

(A+2u)d;.Bk^c^  =  oJ 

where  B  =  |^ . 
a  +  b 

Since  k  =  0,  then  c  =  a'  =  0.   Also  (6.8)  Implies  d'  =  0 

so  that  d  =  A  =  constant. 

o    o 

For  n  ^  1  the  general  solution  of  (6.8)  is 

c     =   {X  +2[i){A     cosh  k  z  +  B^   sinh  k„z) 
n  n  n         n  n 

d„  =  B(A      sinh  k  z +  B„   cosh  k„z ) 
n  n  n         n  n 

where  A  and  B  are  constants. 

Next  Integrating  the  system 

a'  +  k  b  =  c 
n   n  n    n 

b'  +k  a  =  d 
n   n  n    n 

we  obtain  b  =  A  z  +  B  and  for  n  ^  1, 
o    o    o  — 
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^n  =  4  B^z  +  C^)slnh  k^z  +  (^  A^z  +D^)cosh  k^z  , 
A 


n 


2  V-V^^nJ^i-^  V 


+ 


n 


4  B  z  -C  +^^ip5  A 
2  n    n   2k   n 


cosh  k  z 
n 


,  2vLb 


where  A  =  X  +  ^^"^  and  the  C  's  and  D  's  are  additional  constants, 
a  +b         "       ^ 
Next  we  turn  to  the  boundary  conditions-   From  (6.2)  we 

see  that  the  perturbed  boundary  conditions  are 


/ 


(6.9) 


S 


q21  =  q31  ^  0 


•1 

u-^     =   -eR 


> 

i 


for  r  =  R 


i  Q^5  =   0      (1=1,2,3)      forz=±h 


Using  (4.5)  and  (6,5)  the  non-trivial  boundary  conditions 


become 


^^■^  =  2y.  5       f-T^^—  c     -k  b     la     =   0     for  r  =  R   , 
' \S  4-b     n       n  n  J  n 


u       =  aR  + 
.13  ^ 


a  a     =   -eR     for  r  =  R   , 

n  n  * 


Q^  =  y~  (Bc„  -  2uk  b„)a     =0     for  z  =  ±h   , 


n33  _ 


Q^^  =   2Aa+  (A  +  2M.)b  +S         [(A  +  2uL)d     -  2p.k  a    ]3     =0     for  z   =  ± 

^^    '^ —  ■•     ^  n  ^^  n  n  n 


=  ±h 


Next  we  choose  k  so  that  k  R  is  the  n-th  positive  zero  of 
J,,  and  we  observe  from  (6.3)  that  a  =  -e.  Then  the  boundary 
conditions  for  r  =  R  are  trivially  satisfied. 

Observing  from  (6.3)    that   2Aa+(A+2n)b  =   0,    the   remaining 
boundary  conditions  are 


i:-.:^) 


<.     < 


2k 

I  Bc„  -  2|j.k  b  =  0      "l 
n   ^  n  n  I 

(6.10)   /  S  for  z  =±h  ,   n  =  0,1,2,,..  . 

I  (A  +  2n)d^-2nk^a^  =  oj 

For  n  =  0  the  first  equation  of  (6.10)  is  trivially- 
satisfied  and  the  second  becomes  A  =0.  Hence  b  =  0  and 

o  o 

b  =  constant.   The  value  of  b  determines  a  rigid  displacement 
along  the  axis  of  the  cylinder,  so  without  loss  of  generality 
we  may  take  b  =  0. 

For  each  integer  n  ^  1,  (6.10)  gives  us  four  linear 
homogeneous  equations  in  A  ,  B  ,  C  ,  and  D  .   In  order  that  the 
perturbed  buckled  solution  differs  from  the  perturbed  simple 
solution,  the  determinant  of  coefficients  of  this  linear  system 
must  be  zero  for  some  n  *■  1.   This  condition  is 


^    slnh  2k  h  \ 
2(A  +  ^)  (^;i.___^J 


^              slnh  2k  h 
±A-  (5A  +4m.) ^j^ 

for  some  n  ^  1  where  the  proper  sign  is  to  be  chosen.   If  the 
upper  sign  is  chosen,  it  will  be  seen  that  the  thinner  the  plate 
the  greater  the  critical  strain  will  be  which  is  physically  very 
strange.  Also  the  critical  strain  corresponding  to  the  lower 
sign  is  always  less  than  that  corresponding  to  the  upper  sign. 
For  these  reasons  we  expect  the  physically  realistic  choice  to 
be  the  lower  sign,  although  it  Is  interesting  that  two  solutions 
for  e  exist. 

Now  let  e  be  defined  by 


g;.;i»J^ 


A.T 
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p 

/  sinh  2k  h          -^ 

n 

Sinh  2k^h 

A  +  (  5A  +  4n) ^^^ 

n 

(6.11)        e„  =  ...hV  v. .     n  =  1,2,3,...  . 


o 


Then  e  must  be  one  of  the  numbers  e  for  buckling  (i.e.  for  A  , 

n  *^  n* 

B  ,  C  ,  D^  not  all  zero),  and  the  series  in  (6.4)  have  only  one 
n   n   n 

o 

non-zero  term  since  e  can  equal  only  one  of  the  e  (note  that 
e  ^   e  for  n  ^^  m).  Hence  the  n-th  mode  of  buckling  is  given  by 

u""-  =  (l-e)r+  6a  a  +0(5^) 

n  n 


(6.12) 


^^=  (i^lrlr^^  +  ^Vn-^o^^") 


where  the  e  appearing  in  a  and  b  Is  replaced  by  e  and  e  is  an 

unknown  function  of  5  having  the  property  that  e  — >  e  as  6  — >•  0. 

2h 
Let  D  =  -^  (the  undeformed  thickness  to  radius  ratio)  and 

j  tie  the  n-th  positive  zero  of  J-,.   Then  (6.11)  can  be  rewritten 


as 

sinh  Dj 


n.l 


^^•11'^        ^n  =  sinh  Di 

0-+  (2  +  5-)  ^. — ^ 

As  n  — >  oo  ,  e  — > -pr-  so  that  all  the  modes  of  buckling 

start  before  the  deformed  radius  is  as  small  as  -g^  R. 

Let  T  be  the  pressure  per  unit  deformed  area  at  which 
n 

the  n-th  buckled  mode  can  start.   Then  (see  discussion  of 

pressure  on  unbuckled  cylinder) 

Ee 
(6.13)        '^n  =  (l-e^)(l-'^'T+  ^.re^)   ' 


<  ..  -..  .  : 


•,>;.f:t   a  ixsaj-' 


'?  <«s 


jv'    Li" 


r?ilv' 


n 


•••p..-<w 


^ —  .'■ 


26 

As  n  — >  CO  ,   T     — >•  T^    =  ^ ■ — rr  so  that  all  modes  of 

"     °°    {2-0-)  (1+0-)^ 

buckling  start  before  the  finite  pressure  T   is  reached. 

As  D  — >  00  with  n  held  fixed,  e  — >  7=t^  and  T  —>  T  . 

'      n    2+ar  n     oo 

Thus  the  intervals  containing  all  the  e  and  T  values  shrink 

down  to  the  limits  of  the  sequences  [el  and  [T  1  respectively 

as  D  — >  QD  . 

To  compare  our  result  with  that  of  the  thin  plate  theory 

we  observe  that 

p2.2 

^n  =  T?n:TTT^°^^''^^n) 
and 

n  1-r    n     n    ^^(i  -o-^)      " 


The  force  per  unit  deformed  circumference  is 


2hbT^  =  n   +0(h5j;)  =  P(^)  +0(hV) 


2Eh^j„         c:  h        J'n  2      ah 


2Eh^ 


where  F  =  — a  is  the  flexural  rigidity  of  the  plate.  Thus 

5(1-^  ) 
to  lowest  order  the  critical  pressure  here  agrees  with  that  of 

the  thin  plate  theory.   In  some  respects,  however,  the  two 

theories  are  quite  different.   Although  the  critical  loads  agree 

to  lowest  order,  as  n  — >  co  the  thin  plate  critical  loads  — >  oo 

whereas  the  thick  plate  critical  loads  approach  a  finite  limit. 

This  means  that,  no  matter  how  thin  the  plate  is,  the  critical 

loads  for  the  tv;o  theories  differ  by  arbitrarily  large  amounts 

for  large  enough  n,  i.e.  for  a  buckling  mode  of  higher  order. 
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Figure  4  shows  the  results  for  this  problem  graphically 

f 

for  cf  ■=■   .3'   The  critical  strains  e  are  plotted  versus  (Dj  Y 
where  we  recall  that  D  is  the  undeformed  thickness  to  radius 
ratio  and  j  is  the  n-th  positive  zero  of  J,  .  The  graphs  are 
shown  for  both  the  thick  plate  and  thin  plate  theories. 


7.  Problem  B 

Cylindrical  coordinates  are  used  again.   Since  the  axial 

_     "5 

displacement  is  taken  to  be  zero,  we  have  x-,  =  u-^  =  z.   In  the 

unstrained  cylinder  we  have  "^^  —  "^  —  ^9   where  r,  and  r^  are  the 
Inner  and  outer  radii.  The  height  of  the  undeformed  cylinder 
does  not  effect  the  results  and  is  arbitrary. 


>/x. 


-5=.X/ 


Fig.  5 


The  boundary  conditions  are 


(7.1) 


i  T^  =  0   (i=  1,2,3)     for  r 


T^  =  _TN^   (1  =  1,2,3)   for  r  =  r. 


-•Vl 
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where  T  Is  the  hydrostatic  pressure  per  unit  deformed  area 
applied  to  the  surface  r  =>  r^  and  N  g.  is  the  unit  outer  normal 
vector  to  the  deformed  outer  boundary. 

For  r  =  v^,   dS  =  r  dSdz,  N^  =  i,  Ng  =  N  =  0.   Hence 


Also 


T^dS  =  Q^-^'n  dS  =  rQ^-'-dSdE  . 


N^iidS  =  If  K  ll  dSdz  =  u^l^^j  X  u^|^4  dedz 
"  ^ijk^  2  P^^  g^  dSdz  . 


Thus  we  also  have 


and  hence 


T^dS  =  -TN^dS  =  -Te^'^'^P.^P,  ^  dedz  , 

j2  k3 


oil  =  .  T  ijkp  p 

^      r  ^    j2^k3  • 


After  similar  considerations  for  r  =  r.,  the  conditions  (7.1) 

become 

I  Q^l  =  0  (i  =  1,2,5)    for  r  =  r^ 
(7.2)      \ 

L0>^  =  -  \   Te^^^Pj2Pk3  ^°^  ^  =  ^2   • 

We  look  for  a  simple  solution  in  which  the  displacements 
are  entirely  radial.   Hence  the  deformed  cylinder  will  be 
another  hollow  circular  cylinder.   Since  g,  is  directed  along 

the  radius  and  gp  and  g^  are  perpendicular  to  it,  this  means  we 

1  2        ^ 

want  a  simple  solution  of  the  form  u  =  f(r),  u  ■=  0,    \x^  =   z. 

Then  from  (5.3) 


'.  f. 
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/f'    0    0\  'f' 

'^,)  =  /  0    -^f   0  )   and   P(^J)  =1 


Assuming  f  >  0  and  f'  >  0,  then  (p-^*^)  is  symmetric  and 
positive  definite  so  (5.3)  implies  C"^-^  =  g-'"'^.   Also  from  (3.6), 

s,  =  f'+  — f-2.   Substituting  into  (3'9)*  we  obtain  only  one 

I   1   '  b 

non-trivial  equation,  (f  +—  f)   =0.   Hence  f  =  ar+—  for  some 

constants  a  and  b. 

We  now  obtain  s,  =  2(a-l)  and  (from  (3.7')) 

q11  =  2(A-Hi)(a-1)-Bi|  ,    q21  .  q31  .  q  . 

r 

Substituting  into  (7*2)  we  obtain 

2(X-Hi)(a-1)--^  =  -Tra+-^^  , 
r^        rg 

2(A-HJi)(a-l)--^  =  0  . 
^1 

r  2 

'  1  \ 
Let  k  = (  —  )  •   Then  the  solution  to  the  above  equation  is 


2 
1-a  = 


(7.3)       ) 


2(A-Hi)(l-k)+T(l+^  kT 


b  =  2^  (a-l)r?  . 
M,        1 


The  a  and  b  given  by  (7.3)  are  valid  if  T  is  small  enough 
so  that  the  deformed  cylinder  remains  hollow.   The  inside  radius 
of  the  deformed  cylinder  is 


f  -     «-  > 
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f(r,  )  =  ar,  +-S-  = 


(A-ha)2!i  -T)(l--k)r, 


1  r 


I        [|J.+  (A4u)lc]T+ 2^i{>v-Hi.)(l-k) 


From  this  we  see  that  the  inside  strained  radius  goes  monotoni- 
cally  from  r,  to  zero  as  T  goes  from  zero  to  2\i    (notice  that 
0  <  k  <  1).   Hence  we  will  look  for  buckled  solutions  which 
start  at  pressures  T  in  the  range  0  <  T  <  2ii..   For  T  in  this 
range  f  >  0  and  f  >■  0  as  we  assumed  above. 

We  will  look  for  buckled  solutions  which  are  symmetrical 
with  respect  to  a  radius.  Hence  v-je  may  as  well  restrict  our- 
selves to  solutions  symmetrical  with  respect  to  the  x,-axis. 

These  solutions  are  characterized  by  saying  that  u  is  an  even 

2 
function  of  d   and  u  is  an  odd  function  of  9.      Hence  we  write 


(7.^) 


lu^ 


u      = 


b  00 

ar + —  +  5        a   (r)cos  n^ 

^     feu    " 


CD 


\        b    (r)sin  n©        with     b^ 

feu    "^  ° 


=   0 


u-^  =   0   . 


From  (5' 3)   we  obtain 


,11 


=  a  — 7^  +  3        a„   cos  n0 


n 


,12 


p21  ^  r 


-  ^"?  ^n  +  F  ^^^^"^  "^ 
r 

(b'  +-  b    )sin  nd 
n     r     n 


P^^  =  -%  +  -T  +  IZ  (^  a    +  J^  b    )cos  n0 
r       r  r-^  r 

pl3  _  h^3  _  p31  _   p32  _   i33  _ 


=  P' 


=  P>-  =  p^^  =  0   . 


From  (4.2-4)   we  have 


::ry:uora   :.,  ■ 


^'i:^   '-?!:!. •'I 


L  '■■  ■ 
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6^2  =  .c21 


2a 


2a  r 


*  1  i 
(the  other  C  '^'s  are  zero) 


•  ,  1 


s,  =  P  .  =  2a  +  5   (a  +i  a  +  nb  )cos  nS  . 
1     1       ^ —   n  r  n    n 


Let 


(7.5) 


Then 


,12 


r  '1 

i  c  =  a^  +  —  a  +  nb 

In    n  r  n    n 


d  =  b'  +-  b  +-^5  a  . 
n    n  r  n   2  n 


=  -  -^  ZZ  d^  sin  ne 
2a      ^ 

=  2a  +  5    c   cos  n9 
z —  n 


plj 


,2j 


( c  -  nd  )cos  n9 
n    n 


(^r;4^n--%^n)^i"'^^ 


plj 


,2j 


2a 


nd  cos  n9 
n 


4r  5    (d'  +-  d  )sin  n9  . 
2a  ^—   "^  r  n 


Substituting  into  (4.6),  we  see  that  the  variables 
separate  and  the  following  differential  equations  are  obtained. 


(7.6) 


pA+2u)Cj!^-nAd^  =  0 
^(<^n+F^n^  -(A  +  2M.)  -^  c^  =  0 


where  A  =  A  +  2)i  - 


A-Hx 


Since  b  =  0  by  choice,  then  d  =  0  from  (7-5) 
(7.6)  c'  =  0  and  c  =  A  =  constant. 


From 


1     n 


x  -:5 , . 


r.  . 


^'  "'0   -TS 


300( 


B) 


n 


B' 


'     r- 


foJ. 


yrr 


t^     c 


n 


en         ;  -- 


-    Cr,   nlo  (    •'  "^  "^'  n^' 


is^'   «:>rft   JJL 


';   -- 


n     1 


,i:^.r)  .Tto-a:  0  = 
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For  n  =  1,2,...  the  solution  of  (7-6)  is 

where  A  and  B  are  constants. 

Next  we  solve  the  system  (7.5) •   For  n  =  0  we  have 

'1  1  '^o 

a    +—  a     =  A     and  a     =  -o  A  r+ —  where  C     =  constant.     For  n  =  1 
oroo  o2or  o 

the  solution  is 

^^ 


^1  "  B  (3A-A-2ii)A^r^  +  i  (A+A+2m.)B^  log  r+-|  +  D^  , 


B  B  CD 

b-j^  =  ■^[-A+5{A+2iJ.)]A^r  +  |(A-A-2(i)^--|(A+A+2M,)-^  log  r  +  -|-^ 

v;here  C,  and  D,  are  additional  constants. 

For  n  =  2,5,,..  the  solution  to  (7.5)  is 

B  D 

Bn+2A    »      n+1        Bn  +  2A        n      ,    ^   „n-l    .        n 

^n  =   -   U(n+1)  ^n""  "  ^Kn-l)  ^^1  +  ^n"^         +  HS+T  ' 

B(n-t-2)  +2A    .      n       B(n-2)  -  2A  \        „      n.2  _^     ^n 
^n  =       Mn+1)       ^n^ Mn-l)       ^H  "  ^n""         "^  ;n+? 

where  B  =  — |^  and  C  and  D  are  additional  constants. 

a 

°  i  i 
Next  we  turn  to  the  boundary  conditions.   Since  (P  "^ )  is 

diagonal,  from  (7.2)  v;e  obtain 

xil  _   ImpiS^^  P   -iT(e^J^P  P   +e^^^  P  P   ) 

for  r  =  rg. 


Thus  for  r  =  rp. 


•^*"^':(    ^.){i,^^ 


^  r  f       ! 

o  "^        o     t>        o  t.> 

X 

.avln^j-':  .3  9^is   ,-.1  bri&   ^0  vtsdw 

•^  &  ■  ^ 


'    i-j-.:  ■    l-c         .  ni-^ 


.:  0  - 


n^  AS  -  >  "-:       n. 


(1  ' 


'i 


n 


(S4Y)    rr 


'{:>!^S'^^ 


,^1    ^•=  i:tfriT 


3^ 


q11  =   -T(a  +  |)  -  T(a  +  -^)  -T  v^  {nb^ +^  aj^)cos  nO   , 


x21 
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Q  -^  =  -  4?  T  5~~  (na    +  rb    )sin  n0,      and     Q''"^  =  0 
„c       * n         n 


Also   from   (4.5) 


e-^   =    2(A-Hx)4.Bi§  + 


[(A+2n)c^-  2^L(i  a^+ nb^)]cos  n0    , 


q2^  =  XZ  tAd^-f"  ^^n-^F  a^J^sin  ne,      and     Q^^  =  0. 


Since 


2(?v-Hi)(a-l)  -^  =  -T{a  +  -^) 


(from  before    (7.3)),   we   see   that 


2(A-Hx)a-Bi|  ^   .T(^+  b^)  -?.(4+^) 
r  r 


r 


when  r  =  Tp.   Thus  the  boundary  conditions  for  r 


r^  become 


(7.6) 


(  (A+2p.)c^+  (T-  2u.)(i  a^  +  nb^)  =  0  "] 


[Ad^+(T-2^)(^a^+ibJ  =0    _ 
Similarly 


for  r  =  r. 


and  n=  0,1,2,...  . 


(7.7) 


f(A+2M.)c^-  2n(i  a^  +  nb^)  =  0^ 
M^.2^x(-%a^  +  lb^)  =  0    J 


for  r  =  r,  and  n  =  0,1,2,.. 


In  applying  the  boundary  conditions  to  our  solution,  it  is 
convenient  to  express  all  constants  depending  on  T  in  terms  of  a 
single  parameter.  The  most  convenient  such  parameter  seems  to  be 


»  f  r  , 

1 


r ,         ■■1 


.0    -    * 


f!         n         * 


_     fSr 


;>niT  -  -!•  ^o ; 


/,^..b)T-   =  ^-^  (I-i?)(a>f/C)5 


aonxS 


.-,  1 


J.SdJ'    99V:    s- 


<   G^^,^_  ^  ,.,_   ^oys  _^^, 


7C 


;o    ','•'.  .1=1    fJ&XVW 


{  .,on  f 


■;-  1 


\,i.'.ii>ill.:jE.'<;i 


j    0     ■•     I      0;:  '  ^,_;i-/.  3 


!     (T.-^) 


3ffr 


a 
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(7.8)        S  = ^^^^^^  ''     .   .  ^ 

(A-Hi)(l-k)(2ui-  T) 

Then  S  >  0  since  k  <  1  and  T  <  2^..   In  terms  of  thi 
parameter  we  have 


s 


T  = 


2nS(A+ix)(l-k) 
A+2M.  +S(A  +  |x)(l-k} 


o    _  (A-Hi)S+  A+  2|i 
^  ~  (  A  +  2n  )  ( 1  +  S ) 


A  = 


^I(AJf2^i) 

(A-HlJS  +  A+  2n 


R  _  (  A-Hi  )  (  A  +  2|X  )  ( 1+S  ) 
^      lA-Hi)S+ A+  2n    • 

For  n  =  0,  (7.6)  and  (7-7)  give  A^  =  C^  =  0.   Thus  we  have 

a  =  0  in  addition  to  b^  =  0. 
o  o 

For  n  =  1,  (7.6)  and  (7-7)  give  A,  =  B,  =  C,  =  0.   Thus 
D,  XXX 

a,  =  D,  and  b,  = =-.  The  contribution  which  D,  makes  to  the 

position  vector  in  the  deformed  body  is 

D, 
D,  cos  0  g-,  -  ~  sin  9  ^2  "  D^(1,0,0).   Since  this  is  a  rigid 

displacement  along  the  x-,-axis,  without  loss  of  generality  we 

may  choose  D,  =  0  so  that  ^t  =  b^  =  0. 

For  each  integer  n  >  2,  (7.6)  and  (7*7)  are  four  linear 

homogeneous  equations  in  A  ,  B  ,  C^,  D^.   Setting  the 

determinant  of  coefficients  equal  to  zero  and  simplifying  we 

obtain 

[n^(l-k)^+  (l-2k)kA^]S^  +  2[n^(l-k)^  +  k^A^]S+ n^(l-k)^- kA^  =  0 


.1    =;     Ff 


9ff^ 


?*• 


•  D    O      ' 


r\     - 


Sr... 

3^- 


■  i. 


srf;i 


-c  ^Bri: 


iff  -^i 


\   ^r.: 


•  Ci-iCji'u 


•n 


'iXVl!^  V 
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where  X  =  k^+k~^-2.   This  is  the  condition  that  a  and  b  are 
n  n      n 

not  both  identically  zero.   The  parameter  S  was  used  because  of 
the  relatively  simple  form  of  the  above  equation. 
The  solution  of  this  equation  is 


-k^  'h±  kj  (l-k)^+  hh  h„  -  k 


Q n     n  _        n 

S k(l  -^kj+TT 


n        .1^2.  h^  If:  kj(l-k)'^+ith 


n  ^  "^         n 

n^      2 
where  h  =  t—  (1-k)  and  the  proper  sign  is  to  be  chosen. 

Let  r\  =   1-k.  For  a  thin  shell,  n  ~  0  and 
k"=i.m+i%2J.  ,2.n(n-l)(n-2)  ^3^  n(n-l )  (n-2)  (n-3)  ^^  ^oCn^)  , 

k-n=  1+^^  ^Il(n±ll  ^2^  ninilMn+il  ,3^  n(n-H)  (n^2)(n-K3)  ,^0(^5), 
X^=nV^nV^"'^^,^^^NVo(n^)  , 

h^  =  1-n  -  lL^T^2^o(n^)  , 


"n 


h_-k  =  ii-^  ti^+0(ti^)  , 


2 

J(i-k)^  +  4h^  =  2  -n  -^^^-f^  Ti^  +  o(n^)  , 


J(l-k)^  +  4h„  =  2  -  3T1  -  "^  "  ^^  Ti^  +o(ti^)  , 


'n  -  -  -  ^''  -   12 


-k^  -h  +  k]{l-k)^  +4h  =  0(t]^)  . 
n   ^         n 
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Thus  from  the  second  expression  for  S  we  see  that  if  the 


plus  sign  is  chosen  we  will  have  S  =  0(~)  and  S 


•  00  as  Tj  —>   0. 

Since  we  want  S  —>  0  as  n  — >  0  ( so  that  the  critical  pressure 
goes  to  zero  as  the  thickness  of  the  cylinder  goes  to  zero),  we 
see  that  the  minus  sign  is  the  proper  one.  Thus  \ie   define 


(7.9) 


r\ 


n 


n 


k"^  +  k""  -  2 


^  (l.k)2 
n 


{ 


k-  h 


Sn  = 


n 


T 
in 


k^  +  h  +kj  (l-k)^+  4h 
n   -^  n 

2n(A+ji)(l-k)S^ 
A  +  2n+  (A4ixJU-k)a^  J 


for  n  =  2,3,4,. . . 


Prom  (7.9)  it  can  be  shown  that  T  increases  monotonically 

n 

with  n  so  that  there  is  at  most  one  value  of  n  for  which  a^  and 
b  are  not  all  Identically  zero.  Hence  the  n-th  buckled  mode  is 


'  u  =  ar+— +6a  cos  n9  +  0(6  ) 
(7.10)  ^  u^  =  Bbj^  sin  n0  +  O(5^) 


>  for  n  =  2,3,4,. .. 


u'=  z 


\. 


where  T  is  replaced  in  a  and  b  by  T  ,  the  n-th  critical 
pressure* 

In  Pig.  6  we  have  plotted  (l-k)S  versus  k  for  several 
values  of  n.   For  k  =  0  the  slope  is  -1  when  n  =  3,4,5,...  and 
-9  when  n  =  2.   For  k  =  1  the  slope  is  zero  for  all  finite  n. 
From  (7.9)  we  see  that  the  graph  of  T  /2m,  versus  k  is  a  mild 
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distortion  of  Fig.  6,   but  T  /2|a,  depends  on  the  elastic  properties 
of  the  material  whereas  (l-k)S  does  not. 

Tg  agrees  with  Lubkin's  result  [3]  (T^  for  n  >  2  is  not 
derived  in  Lubkin's  work),  and  to  lowest  order 


^n  "  2¥ 


ifi^  (n^.l)(l.k)5+0((l-k)S 


which  agrees  with  the  results  of  thin  plate  methods. 

the  T^'s  are  in  the  range  for  which  the  inside  radius  of  the 

simple  solution  is  positive.   We  notice  that  T   depends  on  the 

relative  dimensions  of  the  hollow  cylinder  which  is  not  the  case 

in  Problem  A.  However  as  r,  — >0  or  r^   — >■  00  (i.e.  as  k  — >  0) 

with  n  fixed,  each  T  converges  to  the  same  limit  o*^  ,  t^  ♦  This 

n       '^  2A  +3(1 

behavior  is  similar  to  the  behavior  in  Problem  A  as  the  thickness 

of  the  plate  goes  to  infinity. 

It  is  interesting  to  consider  the  case  where  the  undeformed 

cylinder  is  solid.   In  order  that  there  is  no  singularity  at  the 

axis  of  the  cylinder,  f  =  ar  and  hence  the  stress  is  constant 

for  the  simple  solution  and  B^  and  D  =  0  in  the"  formulas  for 

n      n 

a  ,  b  ,  c  ,  d  .   The  only  boundary  conditions  are  (7.6),  namely 
that  v;e  have  hydrostatic  pressure  on  the  outside  boundary.   A 
calculation  shows  that  (7.6)  becomes 

(A+2M.)AA^r5  ±  (?.2p.)[[^B  +  ^^^  A]A^r2-(n-l)C^r^-^}=  0  . 


ko 


Hence  A  =  C  =  0  (for  1!  j^  2\x).      Hence  this  procedure  does  not 

uncover  a  buckled  solution  unless  the  hydrostatic  pressure  Is 

2\i,   and  then  all  modes  of  buckling  are  possible,  at  least  if  only 

the  first  order  perturbation  problem  is  considered. 

For  the  hollow  cylinder,  as  the  inside  radius  goes  to  zero, 

T  — »-  -r^t   ,  V^     "^  2ll.  Hence,  as  the  inside  radius  goes  to  zero, 
n    2A  +3ti     ^        '  ^ 

the  critical  loads  for  the  hollow  cylinder  do  not  approach  that 
of  the  solid  cylinder.   Presumably  this  is  due  to  the  fact  that 
in  the  limit  the  hollow  cylinder  has  zero  stress  on  its  axis 
whereas  the  solid  cylinder  does  not. 

8.  Problem  C 

Cylindrical  coordinates  are  used  again.  The  unstrained 


M^: 


(0,R,L) 


z^ 


cylinder  occuppies  the  region 
r  ^  R,  0  <  z  ^  L.   The  condi- 
tion that  the  faces  remain 
plane  and  horizontal  can  be 
written  as  x,  =  0  for  z  =  0, 
x^  =  L(1-e)  for  z  =  L  where  e 
is  a  constant  with  0  <  e  "«  I. 
Since  the  faces  remain 
horizontal  and  since  g-,  and 
go  are  horizontal  vectors, 

the  condition  that  no  shear  stress  is  developed  on  the  faces  can 

1    2 
be  written  as  T  =  T  =  0  for  z  =  0,L.  Thus  the  boundary 


Fig.  7 


conditions  are 


hi 


x^  =  0 


X,  =  L(l-e) 


(8.1) 


3 

t1  = 


2 


for  z  =  0 
for  z  =  L 
for  z  =  0,L 


T^  =  0  (i=  1,2,3)   for  r  =  R  . 


Using  (3«2),  the  boundary  conditions  can  be  rewritten  as 


r  3 


=  0 


u^  =  L(l-e) 


(8.2) 


Q 


13 


q2^  =  o 


for  z  =  0 
for  z  =  L 
for  z  =  0,L 


Q. 


il 


0  (i  =1,2,3)  for  r  =  R  . 


1        2 

We  look  for  a  simple  solution  such  that  u  =  ar,  u  =  0, 


u'  =  bz  where  a  and  b  are  constants.  Then 


/a 


/ 
from  (5.3)  and   (P^^)  =  0 


a 
r 
0 


Since  (P^*^)  is  symmetric  and  positive  definite  (for 


a,b  >■   0),  we  have 


l"-^   =  g'J 


(from  (3.3)), 
(from  (3.6) ),  and 
Q^j  =  [A(2a+b)  -  3?v- Siilg^-^  +2(i.P^J  (from  (3.?'))  . 


s^  =  2a  +  b  -  3 


We  now  easily  verify  that  the  equations  Q  "^  |  .  =  0  are 
satisfied  and  from  (8.2)  v;e  have 
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(8.3) 


b  =  1  -  e  . 


Let  T  be  the  normal  pressure  per  unit  derformed  area 
applied  at  z  =  0,L  to  produce  the  simple  solution.   Then  from 
(3.2)  for  z  =  L, 

T  =  -T^  =  -Q^^  ^  =  E  -4  . 


dS 


? 


Let  P  be  the  total  load  on  an  end.  Then  P  =  ttEeR  . 
As  s  — >  1,  T  - 


E  2 

^  and  P  — >TrER  .  Hence  the  special 


(l+o)' 
strain  energy  density  function  used  here  has  a  consequence  that 

a  finite  pressure  and  a  finite  total  load  would,  as  in  previous 

cases,  collapse  the  cylinder  down  to  zero  height  and  finite 

radius. 

We  consider  buckled  solutions  such  that  the  deformed  axis 

of  the  cylinder  remains  in  the  x-,x, -plane.  We  then  expect  the 

deformed  cylinder  to  be  consistent  with 


^u^ 


(8.4)   {   u^  = 


u5  = 


00   oo 

ar  +  >   ^   a  (r)cos  m9  cos  k  5 
^—pr  '-—TT     mn  n 

m=0  n=0 

OP    00 

\   5   t)   (r)sin  m9   cos  k^z 
■^ — 7^  ^ — 7^     mn  n 

m=0  n=0 

CO    00 

bz  +  "S   5    c  _(r)cos  mQ   sin  k^z 


m=0  n=0 


mn 


n 


where  k  =  H^   and  b   =  c   =0  without  loss  of  generality, 
n    L      on    mo 

From  (5.3) 


^3 


pll 


=  a  +  ^       5   a   cos  m©  cos  k  z 
^ ^^ mn  n 


fl2 
J,21 


m 


(-^  a   +—  b  ^)sin  md   cos  k_z 


n 


a 

T 
r 


2  mn  r  mn 


k  a   cos  m0  sin  k„z 
n  mn  n 


(b'  +-  b   )sln  m9  cos  k^z 
^  mn  r  mn  n 


+  \   >    (-^  b  ^  +  -==5  a„„)cos  m9  cos  k  z 


P^5  =  .  ^ 


r      r 


k  b   sin  mO  sin  k^z 
n  mn  n 


n 


pi 


\    c'   cos  md   sin  k^z 
^ mn  n 


P2 


,33  _ 


^  c„„  sin  m0  sin  k  z 
2  mn  n 

r 


P^-^  =  b  +  >   y""  k„c^„  cos  mO   cos  k  z  . 

J. —  ^ —  n  mn  n 


From  (4.2)  we  have  6^'^=  -C^^  (remember  'b'^^   =  g^^ )  and 
from  (4.3) 


612  ^  J_(pl2.|,21) 
2a 

=  _  i.  3 —  '^    (b'  +~  b   +-2L  a   )sin  m9  cos  k„z 
o  2 —  .;: —  \  "nin  r  mn  ~2  mn  n 

2a  ^ 


cl3  =  __L_(pl3.p31)  =^ 


o    o 

a  +b 


a  + 


S 


5   (c'  +k  a^„)cos  m9   sin  k  z 
'^ —   mn   n  mn  n 


.25  ^  _i_(p23.p2)  =-J-n 

a  +b  a  +b 


(I'  c   -k  b   )sin  mS  sin  k„z 
^2  mn   n  mn  n 


Prom  (4.4) 


s,  =  P^^  =  2a+b  + 


(a'  +—  a   +mb   +k  c   )cos  m©  cos  k  z 
^  mn  r  mn    mn   n  mn  n 


cim  ititi 


^TI^ 


X 


■^    ^'^-<i 


r  5-0, 


5. 


5-.;a  e 


-J.     /' 


kk 


Let 


1   .  1 


V  =  ^mn-^F^mn-^^^^mn+Vmn 


(8.5)       <  ,2m 

Then 

cp-"^   =  -  i  3   5   e   sin  m0  cos  k  z 

^       o  z —  i ^mn  n 

2a 

s^  =  2a  +b  +  ^  m  djj^  cos  me  cos  k^z  . 

From  (5.^) 

plJI  -  "s — "v —  fd'  -me   -k  c'  -k  a   )cos  md   cos  k  z 
^  \.~2 ^ —  ^^mn    mn   n  mn   n  mn'  n 


J 

p20|  _x — ^; —  (e'  +-  e   --2^d   -  k^b   + — g  c„„)sin  m9  cos  k  z 
P   I  -  2 —  2 —  ^  ®mn  r  mn  ^  mn   n  ran  ^  ran  n 

2   ^ 
p3j|  _N — ^ —  fc"  +i  c'  -{^  +  k)c   ]cos  m9  sin  k„z 
P   I  -2 — 2 —  '■^mn  r  mn  ^^   n'  ran'  n 

k 
Alji  _   ^ — r —  rJL  e   + — 2_  (c'  +k  a  „)]cos  ra9  cos  k  z 

C    -  -  2 2 —  L  ^e  ^mn  t  ^  ^   mn   n  mn  '  n 

'  J  2a      a  +  D 

•2J|   r —  r r  ^  If.'     +1  P   )  + 2_(J5.  c   -k  b   )]sin  mO   cos  k„z 

I  j     ^^ ^ 2^     mn     I      iiui       a  +  b  r 

2        ^ 

a3Ji     -     1      ^ — ^; —  fn"    +i  c'     -  (^+ k   )c     +k  d     ]cos  mS   sin  k  z. 
C^''      ='S — r2 2 —  ^Cjj^+^  Cj^     ^"^^'^n^'mn     n  mn'  n 

J  a  +b  ^ 

Substituting  into  (4.6)  we  see  that  the  variables  separate 
and  the  following  ordinary  differential  equations  are  obtained: 
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f  (X  +  2|j.)d'  -  nme   -Bk(c'  +ka   )=0 
mn'^mn    nmn   nmn 


(8.6)  J   ^i(e'  +i  e_)  -  (X+2ti)  4y  d_  -Bk  (k  b 


mn  r  mn 


„<:  mn    n  n  mn  „c:  mn 


m 

—7 


v^ 


B[c"  +i  c'  -  (in^  +  k^)c   ]  -Ak  d 
mn  r  mn  J2       n  mn-"    n  mn 


=  0 


for  m,n  =  0,1,2, .. .  where  A  =  A  +77^  ,  B  =  -^^  . 

a  +b       a  +b 


Eliminating  c   from  the  first  two  equations  of  (8.6) 


we  have 


Hence 


mn  r  mn     2       p,  n  mn 


e   =0  from  (8.5)  since  b   =  0, 
on  on    ' 


e^^  =  (A+2n)(A„  r"^-l+-^) 
mo  mo      „m+i 


for  m  =  1,2,3, . . . , 


e   =  -  A  I  (r*)+  ^  D  K  (r*)   for  m,n  ^   0 
mn   r  mn  m  n   r  mn  m  n        *  '^ 


where  r*  =  J^  k^: 

Bessel  functions  of  the  first  and  second  kinds  respectively 


r  and  I  and  K  are  the  m-th  order  modified 
m     m 


Eliminating  e   from  the  first  two  equations  of  (8.6) 
°  mn  ^ 


we  obtain 


(x+^)(c4<^;,-!j  v-<<' 


mn 


2 

-  Bk  [c"  +i  c'  -  (iB^  +  k^)c   ] 
n*-  mn  r  mn  Z^       n  mn-^ 


=  0 


Adding  this  to  k^  times  the  third  equation  of  (8.6) 
(observe  that  A+B  =  'K+2\i),   we  obtain 


kS 


2 

n 

-r 


d"  +i  d'  -  (IIU+k2)d   =  0  , 

mn  r  mn  VS   n  mn     ' 


and  hence 


<^^^   =  B  ^  +  E   log  r 
oo    oo    oo   ^ 


^mo  -^^^o^""^-^^  ^°^"^=  1'2,3,. 


d   =  B  I  (r  )  +  E  K  (r  )   for  n  ?^  0  ,   m  =  0,1,2, 
mn    mn  m  n    mn  m  n        ^   '       >    >    > 


where  r„  =  k  r. 
n    n 

The  third  equation  of  (8.6)  now  becomes 

2  Ak 

c"     +1  c'     -  (iIL+k^)c        =  -^   [B     I    (7    )  +  E     K    (7    )] 
mn     r     mn       ~2       n     mn         B     ••   mn  m     n         mn  m     n   ■• 

for  n  ?^  0  and 

c"     +i  c'     -4  c        =  0    . 
mo     r     mo       2     mo 

Thus 

n 

for  n  7^  0,    m  =  0,1,2,. . .    . 

From  (8.6)   we  obtain  E^^  =  0,   A^^  =  B^^,   D^  =   -E^^ 

for  m  >  1. 

To  obtain  a  ,  b   when  n  >  1  we  solve  the  first  two 
mn   mn       — 

equations  of  (8.6)  algebraically.  To  obtain  a^^^^  and  b^^  we 
solve  (8.5).   Listing  the  results  we  have 
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r 

1 

°-> 

i     oo 

ss 

oo 

r 

a 


lo 


mo 


2^  A     r^     ^-^^ 


a. 


^lo  1°S  ^-■T  +  ^] 


ran 


>vmH-|ji(m-2)    .        m+1     Am  +  |j.  ( m+2 )     ^mo 
4  ( m+1 )     "mo^        "^         4(m-l)     ^jiwT 

+  a  r*"^""'' + -=^     for  m  >  2 
m  _m+l  — 


Jim  i  T    (r*)A 
p^2  r  -^m^   n'^mn 


n 

A_  ^',- 

n 


*[Ti^^;<r„)-^ri-*rjv^n'iB„„ 


-rorlmf^n'C^-T^^FV'^n'I'm, 


Bk' 


n 


n 

-  ^  K'(r„)F„„     for  n  >  1,   m  >  0 
Bk       m     n     mn  —     '       — 


(8.7)   ( 


on 


'lo 


nra 


=  0     for  n  >•  0  by  choice 


^^'^lo-^^^^^losrA-^ 


,5      r 


X(m+2)  +M.(m+4)    .      ^m  .  A(m-2)  -Hx(m-4) 

TOTT ^mo""   ■*•- — iRmrn 


mo 
„m 


m  _m+2 


for  m  "  2 


1    -r'/_*\«        .  rA    -  ^i,-    X      ?v+2m,   t    /-   \i    m 


mn 


n 


n 


mn 


+   4  -^    Kn,(^n)Fmr,        ^^^    m,  n    >     1 

o  — ^     m    n     mn  — 

n 


n 
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mo 


mn 


=  0  for  m  *■  0  by  choice 


A 


[-#  i:(r^)B_+  I  (F„)C, 


Ek  '•2   m  n  mn   m  n  mn 


n  yrf   {  — 


+  -^  K  (r  )E   +K  (r  )F   ]   for  n  >  1,  m  >  0 
2   m  n  mn   m  n  mn        —  *  "'  — 


oo 


mo 


=  B 


oo 


D 

,m   mo 


d_ 


mn 


^^(\o^  -^)   form>l 


I  (r^)B^^+K  {r„)E^„  for  n  >  1,  m  >  0 
m  n  mn   m  n  mn       —     — 


e   =0  for  n  >  0 
on  — 


In  order  that  we  have  realistic  behavior  at  the  origin 

we  want  the  components  of  u  g, ,  u  gp,  and  u^g,  to  be  bounded 

as  r  — *■  0.  Hence  we  want  a  _,  rb  _,  and  c  _  to  be  bounded  as 

mn    mn      mn 

r  — >  0,  To  accomplish  this  we  choose  the  following  to  be  zero: 

°o>  ^^1'  ^o'  ^m  ^"^  -  2).  \o   ^^  ^  ^^'    (^on'  ^on'  ^on  ^°^  "  -  ^^^ 

^^mn'  ^mn'  ^mn  ^°^  "^  ^  °'  ^^^^' 

VJe  immediately  have  a   =  b  ^  =  c^^  =  0  for  n  *•  1. 
"       on    on    on  — 

Next  we  turn  to  the  boundary  conditions.  First  observe 
that  from  (4.5) 
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a  +b 


O   I 


(ac   -  bk  a   )cos  m9  sin  k  z 
mn    n  mn  n 


Q 


23  _ 


a  +  D  r 


Q-^-'-  =  2{X-Hx)a  +Xb  + 


(2na'  +  M„„)cos  m0  cos  kz 
mn    mn  n 


•  21      


(b 


I    m 

mn  ~ ""?  **mn 


a  _)sin  md   cos  k^z 


n 


Q^-*-  =  '5   5    [Be'  -  (A-A)k  a   ]cos  m0  sin  k„z 
^     ^- —  ^- —  *•  mn       n  mn-^  n 


•        • 


•       • 


and  from  (8.3)  that  b  =  -e  and  2( A-Hi)a +|j,b  =  0. 
Then  (8.2)  gives 


/•  rv.  _  » 


(8.8) 


2iia   +Xd  ^  =  0 
^  mn    mn 


>   ^L-^   a^„  =  0 
mn  „c;  mn 


Be'  -  (A-X)k  a   =  0 
mn       n  mn 


> 


for  r  =  R  . 


From  (8.8)  we  obtain  B^^  =  0,  A^^  =  0  for  n  >  1,  and 

a  =  0  for  m  >  2.   Thus  a^^^  =  b^^  =  0  for  m  =  0,2,3,..., 

Pi 
^lo  "  ^1'  ^"^  ^10  "  "1^-   ^'®  ^^^^  "°^  evaluated  a^^,  b^, 

and  c   when  m  or  n  =  0. 
mn 

For  m,n  >  0  and  r  =  R  equations  (8.8)  give 
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r 


-  5r  [ll(r:)--ir  I„(r:)]A_-4  llir)C 


^n 


m     n 


r "     m     n       mn     13     m     n     mn 
n 


A+2M.X  +  4  (HL  +  F  )]I'(7  ) 


C     rA+2M.     1 
t"^      B      - 


2B 


2 
m 


n  n 

A+2n       A 


B      ■  2B    '  *     _2 '   •  2|i 


+    [(^+1) 
[.  J_  i:(r;)  +  (iIL^+i)I    (r*)]A 


n—   m     n 
2 

7 


(1-IS^)+^]I    (?„)]b^„ 
— ^       2u,     m     n  J    mn 


=   0 


n 


r*     m     n 


Kf     cf     m     n 


mn 


n 


(8.9)    I 


a  —         mn       —       mn       mn 
n  n 


n 


n 


.  A±2!i   [i'(7    )-J-I    (F    )](b 
g—         m     n       —       m     n   M 


mn 


=  0 


a    m 


n 

2Au 


-li  -  -iii-  I    {r*)A        -  ^"^   l'(r    )C 
b  r       ni^    n^^mn       B     "^m^^n^^mn 

n 


+  l(X+^)|i;(r„)-%(|-^r„)y7„)]B^ 


=    0 


n 


We  look  for  a  value  e   of  e  for  which  A   ,  B  ,  and  C 

mn  mn   mn      mn 

are  not  all  zero.   The  determinant  of  coefficients  of  (8.9)  will 

o 

be  zero  for  e  replaced  by  e   ,  and  it  can  be  developed  in  a  power 

mn 

series  in  R.   Also  e   will  be  a  function  of  R   (from  (8.9))  and 

mn 


from  physical  considerations  we  have  e 


mn 


0  as  R 


Then 


^mn  "  0(R^^  ^  =  1  +  0(R^),  b  =  1  +0(R^),  A  =  X  +n  +0(R^), 

B  =  [i  +  OiE^),    and  r  =  k  R  and  r*  =  k  R +0(R^)  (for  r  =  R). 

n    n       n    n 

To  lowest  order  in  R  we  discover  that  the  coefficients  of 
the  first  two  equations  of  (8.9)  differ  only  in  sign.  Hence  we 
replace  the  first  equation  by  the  sum  of  the  first  two  obtaining 
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m+1  t'/^*\   .  rni(m+l)  ,  1 


n 

A  r     m+l  1.(7   )^  [i!i(m+l) +i]i   (7   )}c 
Be       —       mn  — '^  mnj 

n  n 

+  y.[2l^I!l±l+  A     (m(m±ll^-   )]    .(-   ) 
I    ••     B        —       2B         —  n   ■■   m     n 


r  r 

n  n 


'   m  J-  II  A      _ / 


n  n 

for  r  =  R. 

Retaining  the  lowest  order  non-zero  term  in  the  above 
coefficients  we  have 


l^C      .:b^   (m-2)B   +...  =  0. 
2m,  mn   4m,       mn 


Retaining  the  lowest  order  non-zero  term  in  the 
coefficients  of  the  last  two  equations  of  (8.9),  we  have 


m 


m 


>-  mn   p,   mn    Sy,     \i         mn  J 
f -nAjj^  +  [A+2n  -(A-HJ.)m]Bjj^-2(A-Hi)Cj^')  +...  =0 


Since  m  /^  0,  it  follows  that  the  determinant  of 
coefficients  of  (8.9)  is  not  zero  unless  m  =  1.  We  conclude 

c 

that  for  all  e  and  m,n  >  0  we  have  A^^  =  B^^  =  C^^  =  0  unless 

mn    mn    mn 

m  =  1.  Thus  a   =  b   =  c   =0  for  m,n  >  0  unless  m  =  1. 
mn    mn    mn 

For  n  ;^  0,  m  =  1,  and  r  =  R,  equations  (8.9)  become 
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r 


-    1     [l'(r*)  -i  IJr*)]A,     -4  lV(r    )C, 
r?       1     n        r?     1      n   -^    In     B     1     n     In 
n  n 

5-     /  3A+2B    1       A    -   )j.,-   X 


n 


.    ^A+2B  ,    A    .    3A+2B    1  w    ,-    a„  ^ 

n 

(8.10)   J  +A  J_  [    ,,-  ,.J_       ,-  ,, 

B  —     ■•   1     n       —       1     n   ■■   In 
n  n 

.|.  2-|Bi_   [x;(r„,-i  I,,F„,].^  I,(1F„,}b,„  =  0 


n 


n 


-  f  f  Il(r*)A,„.  [^1  Ii(r„)  -I  (i.r„)I,(7„)3B 


b  r 


In 


n 


n 


-B-  lT(r    )Ct      =  0   . 
B      1     n     In 


Replacing  e  by  e..   in  equations  (S.IO)  and  setting  the 
determinant  of  coefficients  equal  to  zero,  we  obtain  an  equation 

o 

for  the  Eiv,'^'  Buckled  solutions  can  occur  only  if  e  equals  one 
of  the  e,  's,  and  the  n-th  buckled  mode  has  the  form 


u  = 


u  = 


u 


5  _ 


ar  +  (P,  +^in  ^°^  k  z)cos  9  ^ 

^1  i 

(-  ^  +  b-|^^  cos  k^z)sin  9    K. 


bz  +c,   cos  Q   sin  k  z 
In  n 


n  =  1,2,3, 


where  p,  is  a  constant.   Since  P,  corresponds  to  a  rigid  dis- 
placement along  the  x,-axis,  we  may  choose  3-,  =  0. 


fij 


--•  "I    '     >;  1 


\A\ 


C 


o 
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From  (8.10)  we  see  that  e.   is  a  function  of  a  and  (^2^)  . 

In  _  L 


In  Pig.  8  we  have  plotted  e  =  e,   versus  (— r-)  for  several 
values  of  cr .   The  graph  for  the  linear  theory  Is  also  shown. 
In  the  linear  theory,  the  critical  value  of  e  is  independent 
of  rf   .  Figure  S  shows  that  this  is  also  nearly  true  in  the  non- 
linear theory  when  ^V-  Is  small.   It  appears  that  e  approaches 
a  finite  limit  as  — t—  — ->  oo  as  was  the  case  in  Problem  A. 


The  numerical  procedure  used  to  determine  the  points  for 

_  p 

Fig.  8  broke  down  for  (— r— )   somevjhat  larger  than  l6  since  the 
procedure  used  for  calculating  I,  and  l'  was  valid  only  for 


small  enough.  However,  there  would  be  no  essential  difficulty 
in  carrying  the  calculations  further  by  other  numerical  methods. 

9.  Problem  D 

Much  of  the  work  carried  out  for  Problem  C  can  be  used 
directly  for  Problem  D,   If  we  let  r,  and  r^   be  the  inner  and 
outer  radii  of  the  undeformed  hollow  cylinder,  the  boundary 
conditions  (8.1)  still  apply  with  the  addition  that  the  last 
boundary  condition  is  valid  for  r  =  r,,  rg.  From  these  we 
obtain  (8.2)  where  again  the  last  condition  is  valid  for 
r  =  r^,  rg. 

The  simple  solution  of  Problem  C  provides  a  simple  solu- 
tion for  Problem  D  with  a  and  b  given  by  (8.3). 

We  look  for  a  perturbed  solution  of  the  form  (8.4).  The 
differential  equations  are  the  same  as  those  of  Problem  C  so 
that  the  general  solution  is  (8.7). 
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The  boundary  conditions  (3.G)  are  valid  for  r  =  r,,  rg. 
From  them  we  obtain  a^^  =  b^^  =  0,  a^^  =  P^,  b^^  =  -P^/r, 
a   =  b   =0  for  m  =  2,5,4,...  with  one  exception.  For  a 
certain  unique  value  of  r^/r^,  a-  and  b,  do  not  have  to  be 
zero.   However  this  solution  might  not  correspond  to  a  buckled 
solution  since  in  the  first  order  perturbation  problem  the 
possibility  of  its  occurrence  depends  only  on  the  geometry  of 
the  undeformed  cylinder  rather  than  on  the  loading.   We  ignore 
this  possible  buckled  solution  rather  than  go  to  higher  order 
perturbation  problems. 

Since  P^  corresponds  to  a  rigid  displacement  along  the 
x,-axis,  we  take  p,  =  0. 

For  m  =  0,1,2,...,  n  =  1,2,3,...,  the  boundary  conditions 

become  linear  equations  in  A^^,  B^,  C^,  D^,  E^,  P^^.   The 

terms  in  A  ,  B  ,  C   are  given  by  equations  (8.9).   The 
mn'   mn'   mn     ° 

coefficients  of  D„„,  E ,  F^„  can  be  obtained  by  replacing  I 

mn   mn   mn  m 

and  its  derivatives  by  K  and  its  derivatives  in  the  coefficients 

m 

of  A   ,  B  ,  and  C   respectively.   The  equations  are  to  be 
mn   mn      mn 

valid  for  r  =  r,  ,  r^-      This  gives  four  equations  in  four  unknowns 
when  m  =  0  and  six  equations  in  six  unknowns  when  m  7^  0.   As 

o 

previously  we  let  e   be  the  value  of  e  which  makes  the 
■^       "  mn 

determinant  of  coefficients  equal  to  zero.   The  e^^'s  are  then 


irnr. 


the  critical  strains  at  which  buckling  can  start. 

From  (8.9)  we  see  that  e„„  is  a  function  of  tr ,  m,   t  > 

mn  -ij 

irnrp 

and  f   .  Figures  9-13  give  an  indication  of  the  dependence  of 

e   on  some  of  these  parameters.   Each  figure  shows  the  graph  of 

"^  r  Trnr^ 

e       versus  100  —  for  (f  =   .3,   one  value  of  — 5 — ,   and 
mn  rp  aj 

m  =  0,1,2,3,4. 
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If,  for  example,  we  consider  a  hollow  cylinder  with 
7rr 
-^  =  1,  then  Figures  9-15  give  e^  for  n  =  1,2,3,^,5, 

m  =  0,1,2,3,^,  and  all  Inner  radii.   In  this  example  we  see  that 

ra  =  1,  n  =  1  is  the  first  buckled  mode  which  can  occur  when  the 

inner  radius  is  small  enough.  Also  as  the  inner  radius  goes  to 

zero,  the  lowest  critical  strain  e,^  approaches  .191  which  is 

the  critical  strain  for  the  solid  cylinder  (see  Fig.  8).  Hence 

the  hollow  cylinder  behaves  like  a  solid  cylinder  when  the  Inside 

radius  is  small.  As  the  inside  radius  becomes  large,  we  see 

that  the  lowest  value  of  e^^  ->  0,  so  the  thinner  the  cylindrical 

shell,  the  sooner  it  v;ill  buckle.   The  graphs  show  that  as  the 

cylindrical  shell  becomes  thinner,  the  value  of  m  for  the  first 

mode  becomes  larger  since  e   is  monotone  decreasing  in  m  for 

mn  ° 

fixed  n  and  r^/r^  =  1,  and  the  value  of  n  for  the  first  mode  also 
becomes  larger  since  the  graphs  become  flatter  at  r^/r^  =  1  as 
n  increases.   Hence,  as  we  would  expect,  the  first  buckled  mode 
is  more  crinkled  as  the  cylindrical  shell  becomes  thinner. 

We  do  not  mean  to  imply  that  the  above  conclusions  hold 
for  all  holloxv  cylinders.  For  example,  if  -j^   =  5,  Fig.  13 
shows  that  e  -j^  <   e,,  for  inner  radii  small  enough  (although 
rather  large).  This  gives  rise  to  the  possibility  that  hollow 
cylinders  of  large  outside  radius  and  small  inside  radius  do  not 
buckle  like  a  solid  cylinder.   Of  course  we  expect  such  cylinders 
to  fail  by  crushing  before  they  buckle  anyhow. 


57 


II 


II 


u 


e 
m 
o 


o 


^ 

^ 


/ 


/ 


/ 


/ 
/ 


/ 


/ 


/ 


t<> 


/ 


o  / 
II  / 
e 


y 


/ 


/ 


/ 


/ 


/ 


/    / 


1      1 

/  / 

i    ; 

! 

i 

;   / 

<      1 

r-\\    OJ 

o 
o 


'^/ 


<^ 


/ 


/ 


o 
o 


ON 


o 

03 


O 


o 


o 


II 


o 


O 
OJ 


O 
i-l 


U 


ON  CO 


vD  ir\ 


KN  OJ 


Fig.   9 


58 


o 

o 


n 
b 

CVJ 

II 


u 
c 


e 
a 
a> 
m 

3: 
o 


o 


(1) 


ON 


CO 


c— 


M3 


ir\ 


K\ 


OJ 


Pig.    10 


59 


cu 
u  \u 

o 

o 


A 


r<> 


cvj 
U 
C 


B 
a 

fQ 
o 


o 


U) 


Fig.    11 


60 


o 
o 


II 
b 

tl 


CM 
U 
C 


a 
n 

o 


o 


u)       O 


o\ 


CO 


Fig.    12 


61 


rH     CM 

o 

o 


A 


II 
in 

I! 
CM 


a 

PQ 
O 


o 


/ 


/ 


in 


o 
o 


// 


^  /  /'/ 


/ 


/  /' 


/     /     / 


'^      !      /  / 

/     /  //■ 

hi  h    Mh/  h 

// 


O 


o 

CO 


o 


o 


o 
in 


o 


o 


<— — 
to       o 


ON 


CO 


t- 


VO 


r<A 


cu 


Fig.    13 


iS-i 


ffi 


<; 


62 


10.  Problem  E 

Spherical  coordinates  are  used  (see  Section  5).  The 
unstrained  region  is  r,  _^  r  _^  rp.   The  boundary  conditions  are: 

I  T-  =  -TN^  for  r  =  rg 
(10.1)         <^ 

j  T"^  =  0     for  r  =  r, 

where  T  is  the  hydrostatic  pressure  per  unit  deformed  area. 

For  r  =  rg,  dS  =  r^  sin  6   dSdcj),  N^  =  1,  Ng  =  N,  =  0.   Thus 

T^dS  =  Q^^.dS  =  r^  sin  9  Q^""-  dSdcj)  . 

J 

Also  _,    ^ 

^'Sl-^S  =  If   ll  dW*  =  "^  I2  Sj  X  u''!  ^  4  d9di 

Thus  v/e  also  have 

T^dS  =  -TN^dS  =  -Te^'^^P-oP,  ,  d0d<|) 
and  hence 

r^  sin  9  ^^  ^^ 

After  similar  considerations  on  r  =  r, ,  the  conditions  (10.1) 


become 


(10.2) 


^''  =  -  :24-n-^  ^'''''^J2\3  ^°^  "  =  "2 


Q"^-*-  =  0  for  r  =  r^^ 


We  look  for  a  simple  solution  having  spherical  symmetry. 
That  is,  a  solution  of  the  form  u"""  =  f(r),  u  =  u^  =  0.   Then 
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(P\)  = 


f '   0 


0   ^f 


0 

0 


f '    0 


(pij) 


0   -f. 
r 


0 


0 


0 
0 
f 


\ 


r^sln  a 


^ifi-^iy   3i  =  f'+|f-3,  p 


'^'  -|:(^'^l^)>  P''l  =P^' 


=  0, 


0    ^^      "        'J        J 
Hence  (3.9)  reduces  to  only  one  non-trivial  equation,  namely 

Thus  f  =  ar  +  — ^  where  a  and  b  are  constants.   It  follows 
r 

that  s,  =  3(a-l),  Q^^  =  oP'^   =  0,  Q^^  =  ( 5^  +  2y.  )(a-l)  --^ . 

The  non-trivial  boundary  conditions  (10.2)  are 
Q^^  =  -T(-)  for  r  =  rg  and  Q   =  0  for  r  =  r^.  These  become 


(10.3) 


(3X+^)(a.l)-iiti|=.(a+^) 


-2 


(3A  +  2n)(a-l)-^  = 


from  which  we  obtain 


r       r^ 
b  =  -j^  (3A+  2u)(a-l) 


(10.4) 


< 


T  = 


(3X+2n)(l-a)(l-k) 
[a  +  ^^^  k(a-l)]'^ 


r  X  3 
where  k  =  (— }  so  that  0  <  k  <  1. 

The  inside  radius  of  the  strained  sphere  is 


f(r^) 


ar  +4 
^1 


■^ 


ar,  +-Tr^  (3A+  2n)(a-l) 
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4ii 
Since  we  want  this  to  be  non-negative,  we  have  1-a  ^  i ,  t^a — r  . 

Since  also  T  ^  0,  we  have  1-a  ^  0  from  (lOJl-).  Therefore  we 

consider  the  range 

(10.5)  0<  l-a<  ^(,V,^^)  ■ 

For  the  above  range  of  a,  T  decreases  monotonically  from 
f  1  kW  ^x  +2u.)   ^°  zero  as  a  increases  monotonically.  Thus  as  the 
hydrostatic  pressure  increases  monotonically  from  zero  to 
( I'^k  W  '^'k  +  Li.  > '    ^^®  inside  strained  radius  decreases  monotonically 
from  r,  to  zero.   Consequently  buckled  solutions  which  start  at 
pressure  T  in  this  range  are  the  only  ones  considered. 

We  will  look  for  buckled  solutions  which  are  rotationally 
symmetric  around  the  x^-axls,  i.e.  for  solutions  of  the  form 
u-^  =  u^(r,e),  u^  =  u^(r,0),  u^  =  0. 


We  write 


(10.6) 


u^  =  f  +  YH   ^n^^^°'n^®^ 


[  u^  =  m  b^(r)Pn(^) 


and  attempt  to  choose  a  and  P  so  that  we  have  separation  of 
variables  in  (4.6).  Then 
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p21  =  P^   =  y~  (b'  +i  b    )3 
1      -^ —      n     r     n     n 

p25  =  0 
Pl=  0 
P2.   0 

P'i.=IZ[(a;;.|a;-^a^)a^ 

p2j|      ^  ^  [(b"+-^  b'  +-%  b    )P 
^^ "■      n     r     n       2     n  *^n 

1.2     n     n       n  g^^2g       ^     n  n 


J 


Ph    =  0  . 


J 

To  have  separation  of  variables  in  the  P  '^  I   themselves, 

,  ,    '  J 

it  is  sufficient  to  have  a" +  a  cot  9   and  p'  +P  cot  9 

n   n  n   n    o 

proportional  to  a     and  also  to  have  P"  +p'  cot  9 ~   and  a' 

^  "   "       sin'^e      "^ 

proportional  to  P  .   It  is  also  to  be  hoped  that  a  choice  of  a 

and  P  which  accomplishes  this  will  lead  to  a  separation  of 
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variables  in  (4.6).   These  proportionality  conditions  are  met  if 
we  choose 


(10.7)  -I 


a^(9)  =  P^(cos  0) 


^  for  n  =  0,1,2, . ., 


i  P  (0)  =  aj^(0)  =  -sin  e   P^(cos  9) 


n 


n 


where  P  (x)  is  the  Legendre  polynomial  defined  by  Pq(x)  =  1, 

P  (x)  =  — i—  ~r  (x^-l)'^  for  n  =  1,2,3, ...  .   Then  u^  will  have 

"      2  ni  dx 
a  simple  zero  for  0  =  0,  tt,  but  this  is  desirable  since  we  expect 

the  solution  to  have  that  property. 

2 
The  differential  equation,  (1-z^)  ^-n  -  2z  ■^+n(n+l)u  =  0, 

dz      ^^ 
for  the  Legendre  polynomials  can  be  written  in  the  form 

a"  +  a'  cot  0  =  -n(n+l)a  .   Since  p  =  a'  from  (10.7),  we  have 
n   n  n         n    n 

p'  =  a"  =  -3  cot  0  -  n(n+l)a  ,  and 
n    n     n  n' 

p"  =  -p'  cot  0  +  p  csc^0  -n(n+l)P  .   We  list  these  as 
n     n         n  n 


(10.8) 


fa     =  p 
n    n 


p'  +  p  cot  0  =  -n(n+l)a 
^n   n  n 

^  a"+a'  cot  0  =  -n(n+l)a 


P 


n 


p;;  +  p;cot0-      ^ 


2 n(n+l)p^ 


These  will  be  used  to  replace  derivatives  of  a  and  6  by 
expressions  involving  a  and  P  themselves. 


For  convenience  we  Introduce 
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To 


=  b    +  —  D r>  a 

n  n     r     n     „2     n 


(10.9) 


^^n  =  <4^n-^("-^l)^n 


where 


Then 


I  e^  =    [h(As^-  2ii)  +  2n]c^ 


h(r)   =  ^r; Y^  = 


f  '  +  ^  f        2ar^  -  "S 


plj 


n  ,  2 


P-I3.  =  Ii:[a;;^£a;-    ^ 


n  +n+2  ^     .  2n(n+l)   ,     , 


[dj!^+n(n+l)c^]a^ 


f,2j 


r>."  J.  ^  I,'      n  +n-2     k    j.  2     ^    ip. 


n     r     n     „2     n     n 


5 —  [c'+-i  c    +Ak  d    ]p 

^1 '•10         v>         »-i  U>         n  ■*  ^ 

r 


C^j    =    -CJ^ 


C^^  =  h(P^^-P^^)   =  -h 


A13  ^  v,^pl3_;31^   ^ 


c  P 
n  n 


^   using  (4.2)   and   (4.5) 


C-"^  =  h(P-'^-  ?''-')  =  0 


0^5 


^  (P 
2f 


23.p32j 


=   0 


J 


o 


=   3(a-l) 


s^  =   3a  +  >        d^a^     using   (4.4) 


a^j 


hZZn(n+l)c^a^  ^ 


A20 


.5j 


(^+|)^h^n^^n       ^   ^^^""S   (5.8) 


C^^l.    = 


0    . 


n'-'^n- 


£f  ^,: 


('  ,01 ; 


f  f 
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Substituting  into  the  differential  equation  (4.6),  we 
obtain  two  non-trivial  eqiiations,  and  we  discover  that  variables 
do  separate  yielding 


(10.10)       y 


(A  +  2|i)dWn(n+l)ej^  =  0   for  n  =  0,1,2,... 


L   r 


^  A^  d  +e'  +-  e  =0    for  n  =  1,2,3,... 
'    n   n  r  n  '  '  ' 


(since  P  =  0,  we  cannot  claim  the  last  equation  for  n  =  0). 
Solving  (10.10)  we  obtain  for  n  =  1,2,3,... 


r 


d  =  constant 


nB.. 


(10.11)     ] 

S 


^n  = 


^n  = 


S  = 


A  r"-l+_n 
n      „n+2 


h(As^-  2g,)  +  2u,  '  ^ 


r 


Next   solving    (10.9)   vje  obtain 


1      ^    ^n-1    r  n     r     n    .^ 

a     =  7s — r=r  \  nr  /     r— ?s dr 

n       2n+l  I  J 


(n+l)c    +—  d 
n     "     - 

3=^ 


._n+l      Hnc     -id    )r"+^drl 
n+2    J  n     r     n  ) 


(10.12)      I 


\-dr.[^''-'f 


Fir^ 


dr 


n+3 


(nc    -id    )r"^^  dr 
n     r     n 


] 


for  n  =   1,2, 3, • • •    • 

Let  H(r)   = 
Then 


(As^  -  2ii)h(r)  +  2|i 


and  L^(r)   =   /     H(p)p"  dp. 
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^n  = 


-^  t(^2n^ 


n-1 


^S)  .i^  + 


n-1 


n+2  2n+l  2(2n-l)  (  A+2^x  )r 


]nB 


+   "   Q  r""^   ^"^^       ^ 


(10.15)  ^ 


^n  = 


2n+l  n 
n+1 


2n+l  ^n+2 


,n-2  ,  n   ^2n+2    (n-t-Dr" 

r 
r 


^2n+l  ^1^    ■*'2n+l^n+3  "  2(2n+3)  ( A+2n  )  ■■% 


+  [..g^   L  _r"-2  •   ^ 


n 


2n+l  L.2n^    "^  2n+l  ^n+3  "  2(  2n-l)  ( A+2^^  )r"^l^^" 
n-2 


C  r 
.  _n . 

+  2n+l  ^ 


D 


n 


(2n+l)r""^^ 


for  n  =  1,2,3,...  where  C^  and  D  are  constants  of  integration. 
Next  we  consider  the  perturbed  boundary  conditions.  From 
(10.2)  for  r  =  rp  we  obtain 


4" 


.ijk 


CO  o      . 


c  o      o 


-^ (T  P.^P,  , +T  P.^P,  ,  +  T  P.^P,  ,) 

r^sin  9  ^^  ^^  ^^  ^^  J^  ^^ 


1 6^23  II 


r^sln  e 


^2^5     r^sin  0  J2^3^^  ^22^k3^ 


Hence 


Q 


11 


Tf^        Tf   A  Tf  ^ 

"7^  "  'p  ^22'  ^3^,„2^   ^33 


r^sin^e 


=  -   (^£5)-^  EZ[|a„-n(n^l)b„]a„ 


Similarly 
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Q-^^-L  =  ^  ^  (a^-  rb^)Pn  ^"^  ^^^  =  0  f or  r  =  rg. 


^21  ^  Tf 

r' 

From  (^.3)  we  have 


q11  =  (3A  +  2^x)4-i^  + 


(Ad^+2,ia^)a^  , 


q2i 


[2ij.(b^!^+i  b^)+  (As^  -2uL)hc^]p^,  and  Q^^  =  0  . 


Prom 


*        2 
(10.3)  we  have  (3A+2|i)a--^=  -(^1^)  for 


^  -  -V— ^;  lux-  r  =  rg. 


il 


Thus  from  the  above  expressions  of  Q   we  have  for  r  =  rp 


o  o 

[^(A+2M.)d^+(^-ai)[|  a^-n(n+l)b^]  =  0 


(lO.lif)  / 


for  n  =  0,1,2, . .. 


Tf 
n   ■  r 


1 

— 7 

r' 


e„-  (4f-  2vl)(-^  a^-i  b^)  =0   for  n  =  1,2,3,...  . 


Similarly  for  r  =  r. 


(10.15) 


n 


r  n 


n- 


s 


'  (A  +2M.)d„  -  2m. [|  a„-  n(n+l)b„]  =0  for  n  =  0,1,2,. 

for  n  =  1,2,3, . 


^n +2.^(7  ^n-F^n)  =0 


For  n  =  1,2,3,...,  equations  (10.14)  and  (IO.I5)  are  four 
linear  homogeneous  equations  in  A  ,  B  ,  C  ,  and  D  .   Setting  the 
determinant  of  coefficients  equal  to  zero,  we  obtain  an  equation 
for  the  critical  values  of  T.   It  is  now  clear  how  the  critical 
values  of  T  can  be  computed  numerically;  however,  we  have  not 
carried  out  such  calculations. 
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Lubkin  [3]  finds  the  lowest  critical  load  for  nearly  solid 
spheres  and  for  thin  spherical  shells,  and  he  observes  that  the 
result  for  thin  shells  agrees  with  that  found  by  the  usual  thin 
shell  procedures.  The  formulation  given  here  can  be  used 
relatively  conveniently  to  calculate  the  critical  loads  for  any 
hollow  sphere  and  any  of  the  modes  of  buckling  considered  here. 
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